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IT OREWAR D 

The purpose o£ this book is two-fold. First, as a 
mathematician, I'd like to teach some relevant mathematics as it 
relates to the modern world of computer graphics. Second, as a 
lover of these "new" high resolution microcomputers, I'd like to 
present some basic principles of computer graphics. I want to be 
able to provide the mathematical background as we go along so that 
you may continue your education in other books more sophisticated 
than this one. 

I'm going to assume at least some minimal background in 
algebra. Trigonometry would be nice. However, I will give enough 
of the basic Trig background to allow you to survive the 
explanations. However, if you really are serious about computer 
graphics, you should get yourself a high-school level trig book at 
your "local" bookstore. The rest of the mathematics you will 
learn traditionally comes from a variety of courses, even Calculus 
and Linear Algebra! (DON'T PANIC— we're not going to need 
Calculus— not everything in a calculus course is calculus!) 

PROGRAM NOTES 

The programs in the book are written in GFA BASIC. I did not 
make full use of the power of this version of BASIC, as my purpose 
was to demonstrate graphics principles, not the language. For the 
most part, medium resolution is what I've assumed. In most cases, 
low resolution would be OK also. Don't be afraid to experiment! 
I've virtually done nothing with color. Again, the purpose of the 
book was not to create paint programs or do arcade(game) graphics. 
Feel free to add some color if you desire to "jazz" up the 
programs. I didn't want to distract from the purpose of this book 
and make "hard" looking programs by putting in fancy stuff. Don't 
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misunderstand me, I love color--I wouldn't dream of not having a 
color computer—I've owned three of them, from a Compucolor, an 
Atari 800, and now my present 1040 ST. Most of all my programs 
that I do for fun involve plenty of graphics and color! 

I've included almost all the programs in the book on a disk. 
You don't know how much I debated over that! I feel there is a 
great value in typing in programs. I felt I've actually got to 
basically understand many programs during the process of copying 
the programs from a magazine article into my computer. The 
temptation is not to closely look at the programs when they're all 
done for you. I hope you will attempt to write your own programs 
and to modify the book's programs—the only way to learn is to be 
an active participant--programming is not a spectator sport. 

MISCELLANEOUS 

I've included answers to almost all the problems in the 
appendix. I've also included in the appendix a couple of other 
program listings related to the book. (They're on the disk also) 

Have fun. 

Brian Hogan 

Instructor at Highline Community College 
Member of S*P*A*C*E computer club 
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chapter one 
POINTS & LINES 

Let's start out by defining a few terms. These won't quite be 
as precisly stated as some of the mathematical definitions that 
will come later on(but I'm sure you'll forgive me). 

PIXEL: The smallest point on the screen that can be illuminated. 

RESOLUTION: The number of pixels that can be put on a screen. 


This is usually given by stating the number of pixels that will 
fit across the screen horizontally then giving the number of 
pixels that vertically go down the screen. 


The standard 

resolution for 

the ATARI 

ST depends, as I'm 

sure 

you know, on 

the resolution 

you choose 

from 

the desktop. 

Low 

resolution is 

320 by 200, 

medium is 

640 

by 400, and 

high 

resolution(monochrome monitor) 

is 640 by 

400. 

Other machines 

may 


have slightly higher or lower resolutions or they may have several 
different graphics modes with varying degrees of resolution. 
(Usually resolution is sacrificed for increased number of colors.) 
Most of the programs in this book will assume you're in Medium 
resolution(640 by 200) on the ATARI ST. 

PLOTTING A POINT 

To plot a point, there must be a coordinate system established 
for the screen. The usual method on the more popular machines is 
to establish the "origin" (0,0) at the upper left corner. (To the 
dismay of the mathematicians ). This is also dependent on the 
language you might be using. GFA BASIC will assume the origin in 
the upper left corner. TRUE BASIC's location of (0,0) is 
completely in the hands of the programmer. 
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A typical screen might be thought of as in the following 
diagram:(if using GFA BASIC) 


0-x-» 639 



To simply plot a single point is not a very illuminating 

task.(pun intended!!) However, in the interest of learning to 

walk before running a marathon, let me describe a very simple 

program. Don't feel you have to type it in however—as I said, 

its not too exciting. 

1 Interactive input to plotting points 
1 GP1 BASIC VERSION 3 

DO 

ISPUT "Enter coordinates of point to plot";x,y 
PLOT x,y 
LOOP 

' Simultaneously press control,shift, and alternate keys 

1 to exit this infinite loop 

END 

I, on purpose, didn't include any safeguards on your input. 
What happens if you enter numbers that are outside of the range of 
the screen? 

For example, try inputting 400,5 or -2,7 or . 

You will find that the ST's GEM operating system is sophisticated 
enough not give you any errors if you try to plot points off the 
screen. This might not be the case if you tried the same thing on 
other machines. Well, as you can see just simply plotting a 

group of points seems rather "pointless"(I'm sorry). Putting 
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these points in some kind of an order makes things happen. For 
example: 

1 Some older in plotting points 
' GFA BASIC 

FOR i=10 TO 100 
PLOT 1,50 
REIT i 

IIPUT *",a$ ! press Return to end progran 
BHD 

As you can probably see, this will give us a line from (10,50) 
to (100,50). By changing the PLOT line, you can get different 
effects. Try changing the PLOT 1,50 to PLOT i,i . Now you will 
get a diagonal line from (10,10) to (100,100). 

pr oblems 

In the program above, change only the PLOT line to do the 
following:(use only PLOT statements that plot a single point.) 

1. Plot a vertical line 

2. Plot a horizontal line that "moves” from (100,50) to (10,50). 

3. Plot a square such that one corner is at (10,10), and the 
corner diagonally across is at (100,100).(you'11 have to "cheat" 
and replace the single PLOT statement with four PLOT statements) 

4. Try these in place of the PLOT i,50 line.: 

(a) PLOT i,9*RND(1) 

(b) PLOT i,9*SQR(i) 

(c) PLOT i,SQR(2025-(i-55)*(i-55)) 

(d) PLOT i,ABS(55-i) 

(e) PLOT i,45-ABS(55-i) 

PLOTTING A LINE 

If you tried some of the suggested lines in problem 4 above, you 
will now realize that plotting points can give you some 
interesting patterns. Of course in such a simple program you are 
limited to what you can do. Now, if you're just interested in 
graphing lines, then using the built-in command LINE (or DRAW TO ) 
is faster and more economical programming-wise. Again, let's look 
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at a simple program: 

1 Line drawing command In Basic 
* GFA BASIC 

DO 

IHPOT 'Enter endpoints of line(ie. 20,30,40,90)..';xl,yl,x2,y2 
LIKE xl,yl,x2,y2 ! DRAW xl,yl TO x2,y2 performs sane action 
LOOP 

' Simultaneously press Control,Shift,Alternate keys to escape 
EHD 

The ST is sophisticated enough that if you give endpoints such 
that if part of the line is off screen, it will not result in a 

fatal error. This would not be true of all machines and/or 

languages. 

The program above is not inherently very interesting because 

there is no built-in "reason" for plotting the lines. Putting the 
lines on the screen in somewhat of a pattern is what becomes 
interesting. 

Examine the program that follows for an example of plotting 

lines. Its just a "for fun" program. 

■ A 'for-fun' program—nested squares 

' GPA BASIC 

LET t=l 

LET a=0 

LET b=45 

LET s=3 

DO 

FOR i=a TO b STEP s 
COLOR t 
PLOT 10+i,10ti 

DRAW TO 100-i,10ti TO 100-i,100-i TO 10+i,100-i TO 10+i,10ti 
' BOX lDH,10ti,lflO-i,100-i could replace 2 lines above 
' The coord, are diagonally opposite corners of rectangle 
PAUSE 5 ! Causes the display to slov up—GFA is fast! 

BEIT i 
LET t=l-t 
LET a=15-a 
LET b=45-b 
LET s=-s 

LOOP ! Be sure you track this by hand to see vhat is happening. 

' Simultaneously press Control,Shift,Alternate keys to escape 
EHD 


If you run this program you will notice nested 
drawn and erased, over and over. The erasing is 


over and over. 


squares being 
accomplished by 
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changing the drawing color to the background color. The variable 
t "toggles" between 1 and 0 . This allows us to change the color 

each time we repeat the loop. Also the values o£ a,b,and s are 
toggled between two values. 

There are almost countless "cute" little patterns that involve 
line drawing. I won't go into more at this time. If you want to 
try your hand at some patterns try some of the problems below. 

pr ob 1 ems 

1. Change the program above so that there are nested triangles. 

2. Make a string art design similar to the following picture. 



3. Draw a simple picture, such as a house with windows and a 
door.. 

4. Draw a small square and make it move across the screen. 


Expanding on these simple programs for a more serious purpose 
will be the aim of future chapters. 
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chapter two 
LINES & VECTORS 

in this chapter. I'll get down to some serious mathematics. 
We'll examine a couple of ways to draw lines, and in the process 
develop some math. The math may not be new to you, it will depend 
on your background. 

It may be you're asking yourself, "What do I need of knowing how 
to draw lines—I have a convenient LINE(or DRAW TO) command.?". 
Well, that may be so in BASIC, but what if you were put in a more 
primitive environment (such as machine language)? Then it might 
be handy to know how to start creating your own routine to draw 
lines. Also, the mathematics we develop will be useful to us to 
explain other aspects of graphics. 

SLOPE OF A LINE 

If you draw a line on a piece of paper, 
you might ask "What characteristics of 
this line distinguish it from some. other 
line that I might draw?". In looking at 
other lines you might conclude that 
certainly the inclination of the line 
makes it distinguishable. A "flat" line 
looks a lot different from a line at some 
angle. 

Then if you drew two lines at the same 
angle, you then would have to conclude 
that the position of the line makes a 
difference. To describe the inclination 
of a line make the following definition: 
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SLOPE OF A LINE = m = dy/dx where 
dy=y=>-y± and dx=x»-x± 

(xi,yi) and (x 2 ,y 3 ) are any two points 
on the line. dy is often refered to as 
the "change in y" and dx is the "change 
in x". The slope is simply a ratio of 
the vertical distance to the horizontal 
distance. 

A FEW SLOPES COMPUTED 

If the line is horizontal, then the dy 
is 0 in value, hence the slope m = 0/dx = 

0. If line goes at a strict 45 degree 
angle, then the values of dy and dx would 
be the same, hence the slope would be 1. 

If the line is vertical then the slope 
is undefined, since the value of dx = 

0. (Recall division by. 0 is undefined) 

Examples: 

(a) Suppose a line passes through (3,5) and (7,24). Then to find 
the slope of that line, first calculate dy and dx. dy=24-5=19 
dx=7-3=4. Then the slope m=dy/dx=19/4 

Note that it would not matter which point we thought of as being 
(xl,yl) and (x2,y2). If we subtracted in the other order, we get 
dy=5-24=-19 and dx=3-7=-4. But m=-19/-4=19/4 still! 

(b) Suppose a line passes through (3,12) and (7,5). dy=5-12=-7 

and dx=7-3=4. Hence m=-7/4. In this case the slope is 

negative . 

If you would take the time to examine 
the difference in the two examples, you 
will note that in (a) the line was 
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increasing upward, whereas in (b) the 
line was decreasing downward. This will 
produce different ’’kinds" of slopes. The 
"uphill" line will always have a pasltlYS. 
slope and the "downhill" line will always 
have a negative slope. What may be 
slightly confusing later in these terms 
"uphill" and "downhill" is when we get to 
the computer. The y-axis of the computer 
is upside down(from top left corner down) 
from the normal mathematics way of 
thinking of the x and y-axis. 

pr ob 1 ems 

1. Suppose a line L is passing through the following pair of 
points. Find the slope of line L. 

(a) (2,5) and (6,12) (b) (25,14) and (44,2) 

(c) (-3,12) and (-5,-10) (d) (5,6) and (12,6) 

(e) (-5,12) and (-5,-3) 

2. If one point on a line L is (5,8) and the slope is 2/5, give 
the coordinates of 3 other points on L. 

3. If one point on a line L is (-4,5) and the slope is -5/3, give 
the coordinates of 3 other points on L. 

4. If two lines were parallel to each other, what would you 
conclude was true about their slopes? 

EQUATION OF A LINE 

In this section will be developed the 
traditional equation of a line that you 
may have learned in conjunction with your 
high school algebra. 
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Suppose we know a point (a,b) is on a 
line, and that its slope is m. Let 
(x,y) be anv other point on that same 
line. Since a straight line could not 
change its slope any place along the 
line, it would have to be true that the 
slope between (a,b) and (x,y) would be m. 

In other words: 

y - b 

- = m 

x - a 

or y-b = m(x-a) [ point-slope form ] 

Anv point on the line we were talking about would have to satisfy 
that equation.( ie. if we "plugged" in the x and y coordinates 
into the equation, we would get a true statement ) 

Examples: 

(a) Suppose a line has a slope of 2 and passes through (3,4). Its 
equation would be y-4=2(x-3). The point (4,6) is also on the 
line since 6-4=2(4-3) is true. Also (0,-2) is on the line since 
-2-4=2(0-3) is true. If you only knew one coordinate of a point 
on the line, you can put it into the equation and solve for the 
other coordinate. 

(b) Suppose a line passes through (-3,5) and has a slope of -3/4. 
It's equation would be y-5=-3/4(x+3). Suppose you knew the x 

coordinate of a point on the line was 3. To find the y 
coordinate, simply solve y-5=-3/4(3+3) for y. You will get 
y=l/2. 
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ANOTHER EQUATION OF A LINE 

Another useful form of an equation is gotten by considering the 
known point to be the y-intercept . Suppose the line has a slope 

of m and a y-intercept of (0,b). Then the equation of the line is 
y-b=m(x-0) which simplfies to 

y = mx + b t slope-intercept form ] 

The slope-intercept form is often handy to get information about a 
line when all we have is an equation in some non-standard form. 

For example, 3x+2y=9 is an equation of a line. If we solve it 
for y, we get y=-3/2x + 9/2. In this form we can see that the 
slope is -3/2 and the line crosses the y-axis at (0,9/2). 

By the way, the equation Ax+By=C is sometimes called a "General 

Form" of an equation of a line. 

TWO SPECIAL CASES: 

(a) The equation of a horizontal line 

through (a,b) is y=b. 

(b) The equation of a vertical line 
through (a,b) is x=a. 

pr ob 1 e ms 

1. Find the equations of the lines through the given point and 

with the given slope: (a) point is (3,6); m=-3 (b) point is 
(-2,5); m=2/3 (c) point is (-2,-4); m=-2/5 

2. Find the equation of the following lines: (a) y-intercept of 

(0,4); slope of -3 (b) horizontal line going through (-4,5) 

(c) vertical line going through (-3,-5) 

3. Determine the slopes from the following equations of lines: 

(a) 3x+y=4 (b) 2x-3y=9 (c) y=3 (d) y-5=3(x-4) (e) x=4 

(f) -2y-3x-10=0 
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4. Find the equation of the line that goes through the two 
points: (a) (3,4) and (-5,9) (b) (-2,5) and (2,6) 
(c) (-4,-2) and (3,-10) 

5. Find the equation of a line that passes through (3,4) and is 
parallel to a line given by the equation 3x+2y=5. 

6. A line that is perpendicular to another line has a slope that 
is the negative reciprocal of the other slope. Find the equation 
of a line that passes through (-4,5) and is perpendicular to the 
line 3x-5y=2. 

GRAPHING A LINE BETWEEN TWO POINTS ON THE COMPUTER 

Now let's use our new-found knowledge to plot a line on a 
computer screen. The program essentially divides itself into two 
major parts: graphing a line whose slope is less than one, and 
graphing a line whose slope is greater than one. To draw a line, 
we will use the point-slope form of a line, solving it for y: 
y=m(x-a)+b. The slope will be calculated from the two given 
points (xi,yi) and (x^yss) . Then the program will use (Xi,y±) as 

the given point (a,b). At first thought, it should just be a 
matter of using a FOR-NEXT loop from xl to x2, calculate the y 
values, then PLOT x,y. 

If we don't have a steep line this 
works fine, but once we get above a slope 
of 1, moving over 1 in the x direction, 
may give a jump of more than 1 in the y 
direction, hence leaving "holes" in the 
line.( See the figure at the side) So to 
handle that case, if the slope is steep, 
then we just solve the equation for x: 
x=l/m(y-b)+a. Then use a FOR-NEXT loop 
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from yl to y2, calculate the x values. 


then PLOT x,y. 

Since plotting a horizontal line or a vertical line doesn't 
require the use of fancy formulas, they are separate routines. 
Finally, we include the special case of plotting a single point if 
both the given points are the same. It seems like a long program, 
but only a part of it is running at any particular time, hence it 
executes fairly rapidly.( As fast as you can expect BASIC to run!) 


*** plot a line *« 

*» LINEPL0T.V1 »« 

*« 6FA BASIC VERSION 3 *** 
variables used 

xl,yl,x2,y2 end points of a line 
dy,dx change in y and x coord, 
a slope of line 
rs reciprocal of slope 
x,y points on line to be plotted 


DO 

INPUT “Input start and end coord.(xl l yi,x2,y2) “jxl,yl,x2,y2 
LET dy=y2-yl 
LET dx=x2-xl 

IF dx=0 AND dy=0 THEN ! a single point 
PLOT x!,yl 

ELSE IF dx=Q ! a vertical line 
FOR y=yl TO y2 STEP SGN(dy) 

PLOT xl,y 
NEXT y 

ELSE IF dy=Q ! a horizontal line 
FOR x=xl TO x2 STEP SGN(dx) 

PLOT x,yi 
NEXT x 

ELSE ! an oblique line 
LET a=dy/dx 
LET ra=dx/dy 

IF ABS(a)<l THEN ! not a steep slope 
FOR x=xl TO x2 STEP SGN(dx) 

LET y=a*(x-xl)+yl 
PLOT x,y 
NEXT x 

ELSE Isteep slope 
FOR y=yl TO y2 STEP SGN(dy) 
x=ra*(y-yl)+xl 
PLOT x,y 
NEXT y 
END IF 
END IF 
LOOP 

’ Siaultaneously press Control,Shi ft,Alternate to escape 
END 
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You may well ask, "how does this fit in with what we've just 
done?”. The answer to that question will be to develop some new 
mathematics and from that develop a different algorithm for 
graphing a line. So be patient as the new tools are defined and 
developed. The general topic of vectors will give rise to many 
other applications in this course, as you will see. 


In this chapter, we will be talking only about "vectors" in the 
plane. The ideas will easily generalize to higher dimensional 
spaces. 


DEFINITION: A vector v from points 

(a,b) to (c,d) is denoted by 
v = <c-a,d-b>. The values of c-a and 
d-b are called the components of v. 
DEFINITION: Two vectors are equal if and 

only if their respective components are 
equal. 

Examples: 

(a) The vector v from (1,7) to (6,10) 
is <6-l,10-7>=<5,3> 

(b) The vector u from (3,1) to (8,4) 
is <8-3,4-l>=<5,3> 

(c) The vectors in examples (a) and (b) 
are equal since their components are 
equal.( ie. u = v ) 


(d) 

The vector from 

(8,4) to 

(3,1) 

is 

<3-8,l-4>=<-5,-3>. 

This 

is a 


different vector than in example (b). 




Geometrically one represents these geometric vectors as "arrows" 
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from the first point to the second point. By drawing a few 
pictures it soon becomes evident that equal vectors are arrows 
that are parallel to each other and have the same length. If a 
vector goes in a different direction or is a different size, then 
the vectors will not have the same components. 

Be sure to examine the "picture" of a vector that 
has a component of 0. These horizontal and 
vertical vectors can be thought of as the building 
blocks of any vector as you will see later. 

In the previous example, u = v = <5,3>. Note 
that a vector from the origin (0,0) to the point 
(5,3) would be equal to those two vectors also. 

Thus a point (x,y) can be thought of as 
representing a vector <x,y> that is drawn from the 
origin to the point. Because of this one to one 
correspondence between points and vectors from the 
origin, you find many books use the same notation 
for a point and a vector. 

I think this is unnecessarily confusing at the beginning, so I 
will use the notations as established in this book. 

Example: Suppose a vector "begins" at (3,4) and has components of 

<-2,4>. What is the coordinates of the "ending" point? 

Solution: Let (x,y) be the unknown ending point. Thus 

<-2,4> = <x-3,y-4>. Therefore -2=x-3 and 4=y-4 (Equal vectors 
have equal components). Solving these, we get x=l and y=8. So 
the ending point is (1,8).## 

DEFINITION: Let u=<ui,u 2 > be a vector 

and let t be any real number(called a 
scalar) then 

tu is the vector <tui,tu 2 >. 
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DEFINITION: If u=<ui,u 2 > and v=<v lf v 2 > 
are vectors, then u+v is the vector 

<Ui+Vi,U 2 +V2>. 

Note: <a,b> = <a,0> + <0,b> by this 
definition. 

Examine the picture to the right to see 
the geometric interpretation of adding 
vectors. 

As you can see if you connect the vectors 
up, beginning to end, the vector you draw 
from the beginning of the first vector to 
the end of the last vector is the vector 
u+v. 

Note that a single vector is just 
geometrically adding up the horizontal 
and vertical vectors made up from its 
components. 




problems 

1. Let u= <2,3> v=<3,-7> w=<-3,4> 

Find (a) u+v (b) 3w (c) v+w (d) -2v (e) 2v+3w 

2. Find the vector from (3,-4) to (-7,-10) 

3. A vector <-3,7> ends at (4,-3). Find the beginning point. 

4. A vector <4,-3> begins at the point (-4,-5). What is the 
ending point? 

5. A vector <5,2> begins at a point (x,y) and ends at the point 
(2y,x). Find the values of x and y and the two points. 

6. Prove that u+u=2u, for any vector u. 
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Sketch <3,4>, 2<3,4>, -2<3,4>, l/2<3 / 4>. Assume all vectors 
start at the same point. I£ u=<3,4>, make some general 
conclusions about tu. 

8. What would the vector <0,0> represent? 

9. Suppose u=<a,b>. What formula would represent the length of 
the vector u? 

EQUATION OF A LINE(AGAIN) 

Although we did not formally prove the fact, the results of the 
previous section(definitions, discussion and problems) should have 
left you with the following theorem. 

THEOREM: tu is a vector that is vector parallel to the vector u, 
and the length of tu is t * length of u. If t is positive, tu is 
a vector in the "same" direction as u and if t is negative, tu is 
in the opposite orientation. 

(You should try to prove this. Make up a definition of the length 
of a vector. How are you going to get a handle on the "direction” 
concept?) 

OK, let's consider getting an equation 
of a line between two points (xi,yi) and 
(xa:,y 2 ). 

Let u = <X2-Xx,y 2 -yi> be the vector 
from (xi,yi) to (x = ,y = ) . Let (x,y) be a 
point on the line segment between the two 
points. Let v = <x-xi,y-yi > be the 
vector from (x±,y±) to (x,y). 

"Obviously" v is a vector parallel to 
u , in the same direction and shorter 
than u. Thus v=tu, 0<t<l. 
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So, _ _ 

v = tu 

<x-Xi,y-y*> = t<xa-xi.,ya-yi> 

<x-xi,y-yi> = <txs-tx t/ tya-ty».> 

Hence, 

x-xi= txss-txi and y-yi= tya-ty* 
x = xi+(x 2 -xi)t and y = y±+(ya:-y± )t 

The above pair of equations, 0<t<l , are called the parametric 
equations of a line from (Xi,y±) to (x»,y a ) . t is called the 
parameter. 

Example: Find the parametric equations for a line from (3,5) to 

(8,15). 

Solution: x=3+(8-3)t and y=5+(15-5)t, 0<t<l. Simplifying we 

get x=3+5t and y=5+10t, for 0<t<l. ## 


problems 

1. what are the parametric equations of the line segment from 

(2,-3) to (7,4)? 

2. What are the parametric equations of the line segment from 

(-3,-7) to (-50,-8)? 

3. What are the parametric equations of a vertical line or a 
horizontal line. 

4. For the parametric equations x=xi+(x 2 -xi)t and 
y=y*+(y*-yi)t, 0<t<l, what point do you get if t=0? if t=l? 

5. For the line segment given by x=7+8t and y=4-3t, 0<t<l, what 

were the two endpoints? 

6. For the same equations in prob. 5, find the intersection point 

of the line segment with the horizontal line y= 2. Do the same 

for y=-2. 
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7. If we allow t to take on other values besides those between 0 
and 1, what points do we get? Generalize your conclusions in a 
statement. 

8. Use parametric equations to prove the midpoint between 
(xi,yi) and (x 2 ,ya) is ( (xi+x 2 )/2, ( yi+ya) /2). 

GRAPHING A PARAMETRIC EQUATION OF A LINE 

I hope you recognize the relative ease in programming a computer 
to graph a line point by point by these set of parametric 
equations. Basically a FOR-NEXT loop using the variable t, t 
going from 0 to 1. Obviously we need to figure out the 
appropriate STEP size so as to get a solid line. To see what this 
step value is, see what value of t is needed to increase the value 
of x and y no more than plus or minus 1 from the starting point 
(xi,yi) (ie. t=0 for this point). 

I’m sure you don't want me to spoil your fun of writing this 
short program. Stop reading this book and write a short program 
that allows you to INPUT the endpoints of the line segment(within 
the boundaries of your screen coordinates), and the program will 
graph a line between the two points. 
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chapter three 
WINDOWS AND VIEWPORTS 


In this short chapter, we will develop one of the main 
mathematical relationships between the "real world" and the 
computer screen. One main purpose of computer graphics is to 
represent real and mathematical objects on a computer. The 
following definitions are usually used in this topic: 

DEFINITION: A closed interval [a,b] = { xe R: a<x<b } 

DEFINITION: A cartesian product A x B. 

A x B = { (a,b): ae A and be B } 

DEFINITION: tm;n] or ([m:n] ) = { xe Integers: m<x<n } 


Example: 

[l,3]x[-l,l]= { (x,y):xe [1,3] and ye [-1,1] } 
(See the picture to the right.) 



Example: [0;639lx[0;199] represents the ST Medium Res screen. 

[0;319]x[0;199] represents the ST Low Res screen. 

The above are sometimes called Screen Domains. 

Example: RxR (often shortened to R = ) represents the plane. 

(Often referred to as the Cartesian Plane.) 
problems 

1. On a piece of graph paper sketch the graph of the following 


sets: 

(a) E-1,3]x[-2,21 (b) [-3,-l]x[-2,0] (c) [-2;3]x[1;4] 

2. What would be the screen domain of the monochrome ST High 
resoulution mode? 

3. What would you guess the definition of AxBxC was? 

DEFINITION: A window is a rectangular subset of RxR. 

DEFINITION: A viewport is a rectangular subset of the Screen 

Domain. 
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It is important you have the following diagram in your mind to 
keep straight these two definitions. They are often confused!! 
BEWARE— even some texts are at times careless in the use of these 
two words. The “confusion in terms comes from the fact that the 
viewport is the computer's re presentation of the window. Hence 
the little rectangles on the computer screen are sometimes also 
called windows by computer companies( Including Atari, Microsoft, 
etc). However, for the mathematics of the situation, we need to 
keep the two concepts separate. Maybe the following diagram below 
will clarify the situation. 



DEFINITION: The function that takes points.from the window to the 
viewport is called a windowing function. . 

Objects in the viewport may look distorted when compared to how 
they look in the window. However it is normally assumed any 
points in the viewport occupy the same relative position as they 
did in the window. For example, a point in the window that was 
located 1/3 of the way between the two horizontal sides and 1/2 of 
the way between the two vertical sides, would occupy the same 
relative position in the viewport. In other words, 
"proportionality" is preserved. 

With the above in mind, look at the diagrams and the resulting 
ratios. I'm assuming the "normal" real-world orientation of the 
smaller y values below the larger values( just the opposite of 
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most screen coordinate systems which assumes the smaller y values 
are toward the top of the screen). I'm assuming the window is 
[A,Blx[C,D] and the viewport is [XS;XL]xtYS;YLl 


RxR 


A-»B 



X-A XV-XS 


B-A XL-XS 


Screen 


XS 


->XL 


YS 


YL 


(xv,yv) 


viewport 


y-C YL-yv 


D-C YL-YS 


Solving these for xv and yv(the screen coordinates) we get: 

x-A y-C 

XV =- (XL-XS )+XS yv = - (YS-YL)+YL 

B-A D-C 

Now as a double check on these formulas let x=A and y=D. This 
point should be in the upper left corner of our viewport (XS,YS). 
Note, indeed, that xv=XS and yv=YS. 

The formulas simplify quite a bit in certain special cases. For 

example, suppose the whole screen is the viewport( XS=YS=0, XL,YL 

are max. screen resolution ). Then we get: 

x-A D-y 

xv -- XL yv = — YL 

B-A D-C 

Be sure you see how these formulas are derived. Don't try to 
memorize these—your brain is worth more than that! 

problems 

1. How would the mapping formulas simplify if you assumed the 
window was [-F,F]x[-G,G] ( F>0 AND G>0 ) and the viewport is the 
whole screen.( Note that our window will put the origin in the 
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middle of the screen.) 

2. What would the general mapping formula be if you assumed that 
the real world had the same "upside-down" orientation that the 
screen has.(ie, larger y coordinates are below the smaller y 
coordinates) 

A WINDOW-VIEWPORT PROGRAM 

Let me present a program that uses the tools we have so far. 
I'm still going to pretend we don't have any line drawing commands 
in our language. Thus, I'll have to use our line drawing routines 
presented in chapter 2 as a subroutine. I won't actually write in 
the routine, you can substitute in either the routine given in the 
book, or the one you wrote using the parametric equations of 
lines. 

The program will just be a simple demo. First, the user will be 
asked for the window, then the user will give the viewport. Lines 
will be drawn on the screen outlining the viewport. Then the user 
will enter the endpoints of line segments(in real world window 
coordinates) and see the line drawn in the correct relative 
position in the viewport. The program then asks for more lines to 
be input. To keep the program simple, no provision has been made 
to end this infinite loop. Certainly feel free to "jazz" this up 
to your heart's content. 
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1 *« Plot a Line in a Hindow 
’ «* HINDOH.Vl 

’ *« GFA BASIC VERSION 3 

' *« variables used—aain progras 

’ a,b,c,d represents window 

' xs, xl, ys, yl represents viewport 

' xw, yw is point in the windowdocal variables) 

' xv,yv is corresponding point in viewport 
1 e,f,g,h endpoints of line in window 
' xl,yl,x2,y2' endpoints of corresponding line in viewport 
' 5,t,i,dx,dy aiscellaneous variables 
PRINT AT(1,20)| 

INPUT "Enter window A,B,C,D ~[A,B3xIC,D] ",a,b,c,d 

INPUT "Enter viewport XS,XL,YS,YL—EXSjXLTxEYSpYL] ",xs,xl,ys,yl 

f 

f Now draw outline of the viewport 
linedrav(xs,ys f xs,yl) 
linedraw(xs,yl 
1inedraw(xl,yl,xl f ys) 

Iinedray(xl f y5 f x5 f y5) 

f 

DO ! now draw user lines in viewport 
PRINT AT<1,20> 

FOR i=l TO 5 
PRINT STRING*(79," ") 

NEXT i 

PRINT AT(1,20); 

PRINT "window = [ 8 ja| a ,“;b;°lxC"jcj","jdi"] B 

INPUT "Enter ends of line (xl,yl);(x2,y2)--xl,yl,x2,y2 ",e,f,g,h 

convert to screencoordinates(e,f) 

LET xl=xv ’ 

LET yl=yv 

c onvert_to_screencoordinates<g,h) - 
LET x2=xv 
LET y2=yv 

Iinedraw(xl,yl,x2,y2) 

LOOP 

’ Below are the subroutines used 

i 

PROCEDURE Iinedraw(xl,yl,x2,y2) ! these xl,yl,x2,y2 are local 
’ vector aethod—assuaes screen coordinates 
LOCAL t,s,dx,dy 
LET dx=x2-xl 
LET dy=y2-yl 

LET s=l/(RAX(ABS(dx),ABS(dy))+l) 

FOR t=0 TO 1 STEP s 

PLOT xl+(dx)*t,yl+(dy)*t ! paraaetric equations of line 
NEXT t 
RETURN 

r 

PROCEDURE convert_to_5creencoor dinates(xw,yw) 

1 returns values of Xv and Yv, screen coordinates 
LET xv=(xw-a)/(b-a)*(xl-xs)+xs 
LET yv=(yw-c)/(d-c)*(ys-yl)+yl 
RETURN 
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OTHER PROGRAMMING SUGGESTIONS: You might do one or more of the 

following: 

(a) Rewrite the input of the previous program so that one only 
enters , say,A and B, for a window of [-A,A]x[-B,B‘]. Simplify the 
corresponding window function. 

(b) Rewrite program so the entire screen is the viewport. 
Simplify the corresponding window function. 

(c) Allow entry of any endpoints( "real-world" coordinates ) 
of a line, but only the points that would be in the window(or 
viewport) would be graphed. 
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chapter four: 

FUNCTIONS AND 
CRASHING FUNCTIONS 

FUNCTIONS 

For most purposes we can think of a function as a special type 
of an equation that is written in such a way that there is only 
one output for each input. Traditionally, when working with 
numbers, X is often the variable used for the input, Y the output. 
The slope-intercept form of a line, y=mx+b, is in that form—input 
a value for x, and out comes a value for the y coordinate. The 
more general point-slope equation of a line needs to be solved for 
y, in order to think of it as a function: y-b=m(x-a) is 
equivalent to y=m(x-a)+b. 

Mathematicians often think of functions in a more general way 
than just described. The only requirement is that a particular 
"input”(might not be just a number) results in a single 
"output"(again, it might not be just a number). A couple of 
examples that we have met in previous discussions might illuminate 
that remark. 

Example: The parametric equations of a line x=xi+(x 2 -xi)t 
y=yi+(yz-Xi )t can be thought of as function that accepts a single 
input for t, then outputs a single point (x,y). The input is a 
number, the output is a point. 

Example: The windowing function accepts a point (x,y) in the 
window as input, and outputs a single point (xv,yv) in the 
viewport. 

If you're willing to carefully define by what you mean as a 
"single input" and what is meant by your corresponding output, 
you'd be surprised how many things can be looked upon as 
functions! 

Often times it is convenient to employ special function notation 
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in describing your function. One nice feature would be to be able 
to "name" your function. Saying "y=" can be awful 
undescriptive, especially if you're talking about several 
functions at one time. 

The special function notation goes in the following format: 
name(input)=output 

Example: f(x)=2x-5 f is the name of this function, a value for 
x is to be input, the value of the expression 2x-5 will be the 
corresponding output. 

Example: g(t)=(t-5)/(t-4) g is the name of this function, a 
value of t is to be input, the value of the expression (t-5)/(t-4) 
is the output 

Example: p(t)=(3+5t,2-4t), 0<t<l p is the name of the function 
that accepts values of t for input and outputs a point on the line 
segment between the points (3,2) and (8,-2). (Did you recognize 
the parametric equations of x=3+5t and y=2-4t?) 

Now, if you knew the functions that had been defined, and 
someone mentioned function g, for example, you'd know - exactly 
which function being referred to. 

The names f,g,h,etc are often used as kind of temporary names. 
You would not be expected to know what function f was unless you 
had seen it recently defined. Other names are more permanent and 
universal in their usage. 

Some common names are: cos(x), sin(x), log(x), ln(x), .... These 
are names for some trigonometry functions and logarithm functions 
found on most scientific calculators. Sometimes special function 
notation is used, such as |x| is a common symbol for the 
absolute value function. In BASIC we use ABS(x). Another 
special notation is the square, root function, vx . In BASIC we 
use SQR(x). 
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The function notation has another big advantage: namely in being 
able to indicate that a particular input is to be actually used in 
the function. 

Example: Suppose_ I define f(x)=2x+3. Then if I write f(3), 

this is to indicate the value of the function f, when 3 is used 
for the input. Thus f(3) =2(3)+3=9. Similarly, f(l)=5; f(-3)=-3; 
f(-10)=-17; f(a)=2a+3; etc. 

Example: Suppose f((x,y))=(5x+y,2x-y). Then if I write f((3,4)), 
this is to indicate the output when the input is (3,4). In this 
case, f((3,4))=(5*3+4,2*3-4)=(19,2). Thus the input of the point 
(3,4) into the function f is the point (19,2). 

pr oblems 

1. Look in your BASIC manual. What other functions are 
predefined and are part of the language? 

2. What function in BASIC accepts a number and outputs a 
character? 

3. What function in BASIC accepts an input of a "string" and 
outputs the number of characters in that string? 

4. What function in BASIC accepts a pair of numbers for input and 
outputs a lighted "dot"(pixel) on the monitor screen? 

5. Describe how a SUBROUTINE might be viewed as a function? 

6. Let f(x)=3/x. Find (a) f(2) (b) f(-3) (c) f(l/3) 

7. Let f(t)=t A 2+3. Find (a) f(3) (b) f(-3) (c) f(1.2) 

8. Let f(v)=(v,v). Find (a) f(2) (b) f(4) (c) f(-3)) 

9. Let g(c)=(2c,4-c). Find (a) g(3) (b) g(-3) (c) g(p) 

10. Let w((x,y))=(3x-2y,xy). Find (a) w(2,l) (b) w(-3,l) 

(c) w(0.5,0.7 ) 

11. Let f(x)=SQR(x). Find (a) f(3) (b) f(-4) (c) f(16) 

(d) f(6.7) 
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12. Let f(x)=SIN(x). Find (a) £(0) (b) f(1.5) 13. Make up 

a function that gives you the points on a line between (3,4) and ^ 
(-4,5). 

14. Make up a function that takes in a number as input and 
outputs twice that number increased by 7. 

15. Make up a function that takes in a point as input and outputs 
the distance that point is from the origin. 

16. Make up a function that takes in a pair of points as input, 
and outputs the vector from the first point to the second point. 

GRAPHING A FUNCTION 

Graphing a function "usually” implies a function that has 
numerical input and output. The input and output are often 
displayed as a point (input,output) (WARNING:. NOT ALWAYS J!) 

In graphing a function that has a 
single number input and a single number 
output, often the Cartesian(rectangular) 
coordinate system is utilized. The input 
is the x coordinate of the point, the 
output is the y coordinate of the point. 

See the example of graphing f(x)=2x+3 
to the side.(This is obviously a line, 
since it is the form of y=mx+b.) 

In graphing a function like our parametric equations of a line, 
often the input is supressed, and all the graph displays is the 
output(which is a point). 

Sometimes additional labels on the 
graph might indicate a few typical values 
of the parameter t. See the example to 
the side. 
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Later in the book, we will discuss different coordinate systems, 
discuss the concept of parametric functions in more detail, etc. 

For now, let’s concentrate on functions that are of the form 
f(number)=number. Drawing graphs of these functions are often 
rewarding and exciting. We've already looked at the function that 
gives a line. What about functions like y=x a (f(x)=x a ); 
y=sin(x) (f(x)=SIN(x)); y= (x-l)/(x +1); or ..... 

Also, drawing these functions by hand can be very time-consuming 
and frustrating. It involves either knowing beforehand what the 
equation will give you as a result, or else computing a lot of 
points to put on the graph. 

Example: Graph the function y=x* . If 
you don't know what it looks like, you 
need to compute some points. One way to 
organize your work is to make a chart: 

-4 -3-2-10124 


16 


16 



The numbers in the chart give you the points (-4,16), (-3,9), 

.(4,16). Of course this is only a small insignificant 

fraction of the possible points. Hopefully it gives the person 
trying to graph the function a good idea of how ,the function 
looks. If the results are still too rough, pick more points; 
maybe pick some fractional values for x, etc. 

Ok, how could the computer help us? One way is to use it to 
construct our points to be plotted: 

FOR x=-10 TO 10 STEP .25 
y=x A 2 

PRINT "C 0 !x; B , B ;y;") a 

NEXT x 
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(In GFA BASIC you could also do the following:) 

DEFFN square(x)=x A 2 
FOR x=-10 TO 10 STEP .25 
PRINT Tix; ViFN square(x) 

NEXT x 


Presto, in an instant you've got more points than you know what 
to do with!! 

But, obviously, this isn't the end of the story(Did you notice the 
title of the book?). Since we have the commands available in 
BASIC to plot points and we know how to put those points in any 
viewport on the screen we wish, why stop with the generation of 
the points? Let's PLOT those little "buzzards” and have the 
graphing done for us!! 

A GRAPHING PROGRAM 

I'll outline a way to do a graphing program. The program L 
wrote is in the appendix. I don't want to tempt you to take the 
easy way out and not do the program yourself. It's really quite a 
short program. If you think you can do one by., yourself, quit 
reading and DO IT!! If not, read on to see how I did it (I lay no 
claim that it's the best way!). The program I'll describe is not 
really the one I'd like to do. These newer BASIC's do have some 
drawbacks compared to some of the older BASIC's with line numbers. 
For example, I used to be able to say LIST 600, and it would 
display what was in line 600 of my program. I could even do this 
while my program was running. Thus if my line 600 was 
600 Y=COS(X) it would list that trig function. If that was the 
function I was graphing (or about to graph) I could display it on 
the screen while the program was running. I could ask the user if 
that was the function they wanted to graph. If not, I could stop 
the program and have them enter in a new line 600 with a different 
function. Then I could start up the program again, by saying RUN 
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100(for example), and the program would continue at line 100 . I 
could work it in such a way that previous graphs were not 
automatically erased. In other words, I had a kind of an 
interactive function input in my program. I find it almost 
impossible to duplicate this in any simple fashion in GFA BASIC. 
This is unfortunate, in that it is often instructive to compare 
different graphs overlaid on top of each other. Well, enough 
crying for what I "ain't got"!! We'll write a program where you 
must edit the program before you run it in order to get the 
function you want to graph. The program will allow you to see 
that function plotted in different viewports with different 
windows so you can see different parts of the graph on the same 
screen. 

1. At the beginning of my program, I have these statements that 
contains my function. 

LET f$="Cos(X)" 

DEFFN f(X)=COS(X) 

This can contain any function. . I put the function in a string 
also, so I can let user know what function is in the program. If 
you wish to graph another function, stop the program, edit these 
two lines to contain the desired function, then restart the 
program. 

2. Now, back to the design of the program. It might be nice to 
reserve the bottom 4 or 5 lines of the screen for a "text" window. 
We're going to need some place to display the current values of 
the viewport and window, and then ask for input about new values 
of the viewport and window variables. You might consider the 
screen looking, as follows, when the program first starts 
running.( Of course this all just personal taste—feel free to be 
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3. I'd display the current value of the viewport, ask if that was 

ok; if not, then ask for new values of the viewport. Do the same 

for the window. Then for jazziness, draw a box on the screen that 

represents the window.(Use the built-in line drawing commands.) 
Just a simple 

BOX xs,ys,xl,yl command will do the job. (If you look 
this command up in your manual, you'll find you need to give the 
diagonal coordinates of the rectangle.) 

4. Next draw in the x-axis and y-axis if they belong in the 

window(vlewport). This is simply done by determining if 0 is an 
element of [A,B] and [C,D], If the axis do belong, you'll have to 
convert the window value of 0 into its appropriate viewport 
coordinate by use of the windowing function. For example, if 0 is 

between A and B, then there should be a y-axis showing. If you 

give the command 

DRAW -a/(b-a)*(xl-xs)+xs,ys TO -a/(b-a)*(xl-xs)+xs,yl 
The complicated expression is simply what you get when 0 is 

substituted for x in the windowing function. Now do a similar 

thing to plot the x-axis(if there is one that should be showing). 

5. Now graph the function. This is a simple 

FOR x = a TO b STEP ?? !Put in appropriate step value 
LET y=FN fix) iget y-coordinate from our equation 

» CHECK TO SEE IF Y VALUE IN UINDOU 
' IF IN WINDOW, CONVERT TO VIEWPORT Put in appropriate 
» COORD. BY WINDOWING FUNCTION codel! 

PLOT xv,yv 
NEXT x 

That's it! Of course, after this you'll ask if they want more 
graphs. If things are done right, you should be able to pick a 
different viewport and graph the function in a different place on 
the screen. If the user picks fairly small viewports, she could 
have 4 or 5 different graphs on the same screen. It really looks 
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quite nice. 

Well I hope that gets you going. If you want to compare yours 
with mine, check my program in the appendix.(Or examine it on the 
disk) I'm sure yours will be much better! I'll also include on 
the disk a little more jazzed up version that uses the mouse to 
define the position of the viewport. 

Just a closing remark—it would really be nice if put text on 
your axis to indicate the beginning and ending values on your x 
and y axis(ie. window values). The nice thing about the ST is 
that it is easy to mix text and graphics. See your DEFTEXT and 
TEXT command in your manual. 

I decided not to go into detail about that kind of stuff as it 
kind of covers up the purpose of the program and this book. By 
adding all of the "goodies" you could make this program a very 
powerful tool.(Some of the goodies would be interactive function 
input, labeling the axis, color, perhaps a little sound, a screen 
dump to a printer, etc. 
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chapter five 
ESSENTIAL TRIGONOMETRY 




I fervently hope you’ve had some exposure to trigonometry. It's 
going to make this a lot easier. If you don't, I hope what I have 
to say in this chapter will get you to your nearest bookstore to 
buy a trig. book.(Or enroll in a course at your local "college") 
I'll try to present the most relevant parts of trigonometry as it 
will relate to the rest of this book. I only hope I can convey 
the importance of these functions in such a short space. 

EQUATION OF A CIRCLE: 


I'm sure you will agree that a 
reasonable definition of a circle of 
radius "r" and a center at some point . 0, 
is that it is a set of points in the 
plane that are all at a distance r from 
the point 0.(See picture to the side) 

Let's give the point 0 coordinates 
(h,k) and pick some arbitrary point .oh- 
the circle and call the coordinates 
(x,y). By what I said above, the 

distance between these two points is r. 
Using the Pythagorean theorem you get 
the following equation: 

< — h ) = -I- < y — k ) =* = = 



(the difference in the x coordinates give the horizontal distance 
between the two points, similarly the difference in the y 
coordinates give the vertical distance) 

Example: The equation of a circle of radius 5, centered at 
(3,-4), is (x-3)= + (y-(-4)) 2 = 5= or (x-3)= + (y+4)* = 25 
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Example: l£ the equation (x-2) 2 + (y-1) 2 = 9 

is given, you know it is a circle centered at 
(2,1) with a radius of 3. 

A special form of this equation is when the center is (0,0). 
Then the equation is simply 

= nr ^ 



(a) center at (5,5); radius is 10 

(b) center at (-3,6); radius is 2 

(c) center-at (-4,-6); radius is 1/2 


(d) center at origin; radius of 1 

(e) center at origin; radius of 5 

2. Give the center and radius of the following circles: 

(a) x 2 +.(y-2) 2 =36 ' ’ 

(b) (x-4) 2 + (y+3) 2 =17 

(c) x 2 + y 2 = 8 

(d) (x+3) 2 + (y-4) 2 = 64 

3. If (3,5) and (7,9) are endpoints of the diameter of a circle, 
what is its equation? 

4. If x 2 + y 2 =4 is the equation of a circle, and if the 
x-coordinate of a point on the circle is 1, what is the 
y-coordinate?; if y-coordinate is 1.5, what is the x-coordinate? 

5. If (x-2) 2 + (y+3) 2 =10 is a circle, and the x-coordinate of 
point on the circle is 2, what is the y-coordinate? 
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TRIGONOMETRIC FUNCTIONS 

There are two main ways to define the trigonometric functions. 
I'll present both for comparison, as each has their own special 
advantages and application. The two methods are totally 
equivalent though. 


The Unit Circle Definition: 

Let x 2 + y 2 =1 be the equation of a 

circle at the origin, with radius of 
1.(Often called the unit circle). Let t 
be a number representing distance 
traveled along the unit circle from the 
point (1,0).(see diagram to the side) If 
t>0, move counter-clockwise; if t<0, move 
clockwise. The point on the unit circle 
you arrive at has some coordinates (x,y). 
DEFINITION: cos ( t > — >c 
s i ri ( t ) = y 

These are called the cosine and Sine 
functions respectively. 

Example: If t=0 then the point is (1,0). 

Thus the cos(0)=l and sin(0)=0. If 
you've got a calculator, check it out—Be 
sure you're not in DEGree, but in RADian 
mode. Or use your computer. Type in 
? COS(0),SIN(0) 

Example: If t=PI/2 (PI=3.141592654..), 

then you're at the top of the circle at 
(0,1). ( The circumference of the unit 

circle is 2*PI*1=2*PI. The top is 1/4 of 
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the way around—1/4 of 2*PI is PI/2 ) 

So, cos(PI/2)=0 and sin(PI/2)=l 
Example: if t= -PI/2, then you're at the 

bottom of the circle at (0,-1). Thus 
cos(-PI/2)=0 and sin(-PI/2)=-l. 

If you want to see some typical values of sine and cosine 

functions, write a short program: 

PRINT "TVCOSITIVSIHIT)' 

FOR t=0 TO 2*PI STEP .1 
PRINT t,COS(t),SIN(t) 

NEXT t 

To see that these values are indeed on a circle, RUN this 

following program: 

FOR t=0 TO 2*PI STEP .1 
LET x=C0S(t) 

LET y=SIH(t) 

LET xv=50*x*50 
LET yv=-50*y+50 
PLOT xv,yv 
NEXT t 

To understand the computation of the variables xv and yv you 
must realize I'm using a window of [-1,1]x[-1,1]_and a viewport of 
[0;100]x[0;100]. Substitute these'into the windowing function and 
you will get the corresponding equations for xv and yv. 

RADIAN MEASURE OF AN ANGLE 

The numbers t are technically 
measurements along the circumference of 
the unit circle. However, it is common 
to use these same numbers as the measure 
of the central angle.(See figure to side 
of page) This type of angle measurement 
is called RADIAN measure(often 

abbreviated RAD). Thus a RAD measure of 
PI/2_ is equivalent to 90 degrees, PI (half 
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way around the circle) is equivalent to 
180 degrees, etc. I think you see why we 
have to tell a calculator which method of 
angle measurement is being used. If you 
have a calculator^ put it in RAD mode, 
then find cos(PI/4). Now, put it in DEG 
mode, and find cos(45). The answers 
should be the same since PI/4 is 
equivalent to 45 degrees. 

GFA BASIC only allows the Radian method of measuring angles for 
it*s trigonometric functions. This is not true for all versions 
of BASIC. 



SOME TRIG. PACTS AND RELATIONSHIPS 

Since cos(t) and sin(t) are points on the unit circle, it is 
"obvious" that 

-1< cos(t) <1 and -1< sin(t) <1 Also since x 3 + y 3 = 1 

for any point on the unit circle, you get 

(cos(t)) 3 + (sin(t)) 3 =1 or more usually as 

cos 3 ( t ) + s i n 31 ( t ) = 1 

This last relationship is often referred to as the 

fundamental trig. identity 

Example: If sin(A)=3/4 then the cos(A) can be computed as follows 

cos 3 (A)+sin 3 (A)=1 
cos 3 (A)+(3/4) 3 = 1 
cos 3 (A)=1- 9/16= 7/16 
cos (A) =± s/ 7/4 . 

Whether you pick the plus or minus sign for your answer depends 
on which quadrant of the circle A is in. 

The following relations come from the symmetry of the circle: 
cos ( —t > = cos (t ) and 
s±n< —1)= —sin(t) 
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AN ALTERNATE DEFINITION OF TRIG FUNCTIONS 

The following definition is often referred to as the "right 
triangle" definition. Picture a right triangle as below: 

B 
a 

A bC 

DEFINITION COS(A)=t>/c £> ± n ( A ) = <a / c 

Sometimes these are memorized in the following way: 
cos(angle)=adjacent/hypotenuse and 
sin(angle)=opposite/hypotenuse 
This allows for more flexiblity since 
many times the triangle isn't oriented 
exactly like the picture and different 
labels might be used for the sides. 

Note that in the right triangle definition , it is easy to see 
that sin(A)=cos(B) and cos(A)=sin(B) or since B=90-A, 
sin(A)=cos(90-A) and cos(A)=sin(90-A) (we are thinking of the 
angle measured in degrees here. If We were using RAD measure, the 
same identities would be cos(A)=sin(PI/2-A), etc.) 

Note that since the sides are always less than the hypotenuse, 
that the sine and cosine are in general less than l(as we noted in 
circular definition). If you look at the picture to the right you 
should be able to see the reason why the 
two definitions basically give you the 
same numbers—just keep in mind that 
similar triangles have the same ratio for 
corresponding sides. 

You can get the fundamental trig identity from this definition 
also. 
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cos 2 (A)+sin 3 (A) = (b/c) 3 + (a/c) 3 = 
b 3 a 2 a 3 +b 3 c 3 

C 3 C 3 C 3 C 3 

problems 

1. Look at the diagram to the side of this problem. Answer the 

C 

following questions: 

(a) sin(A)= ? (b) cos(B)= ? (c) sin(B)= ? 

(d) sin 3 (A)+ cos 3 (A)= ? 8 

2. Use your calculator(or computer) to find the values of the 
following: 

(a) cos(30 degrees)=? (b) sin(70 degrees)=? 

(c) sin(2.45)=?(use RAD mode) (d) cos(-3.45)=? 

(e) cos(150 degrees)=? 

3. Suppose a right triangle has the hypotenuse of 20 feet and one 
of the angles is 35 degrees. What is the length of the leg 
opposite the 35 degree angle? 




SOME OTHER DEFINITIONS AND RELATIONSHIPS 

Now that you know what we're talking about, let's delve into a 
few more things. What we're doing now won't really depend on 
which definition we use for the basic two trig functions. 

There are four(4) more trig functions that can be defined: 
tan ( t ) = sin(t)/cos (t) ( tangent function) 

cot ( t ) = l/tan(t) ( cotangent function) 

s & c C t ) = 1/cos(t) ( secant function) 

esc C t ) = l/sin(t) ( cosecant function) 

There are two more fundamental trig identities involving these 
new functions. 

1 + tan 3 (t) = sec 3 (t) 

1 + cot 3 (t) = esc 3 (t) 
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(Can you establish these from our definitions and the original 
fundamental trig identity??) 

There are literally hundreds of relationships between these . six 
trig identities. Pick up a trig book to see a few of these. Most 
are not to be memorized, but are just drill type problems to 
increase your familiarity in using the basic definitions and 
fundamental identities. 

A couple of identities that we will need later are these: 

** cos(x+y)=cos(x)cos(y)-sin(x)sin(y) ** 

** sin(x+y)=sin(x)cos(y)+cos(x)sin(y) ** 

These are relatively hard to derive. Drag out that book I told 
to to get!! Actually, only one of those is hard to get, the other 
comes relatively easy from the first one. 

A very useful identity that involves non-right triangles is the 
following: 



a 


= si 23 * -+- b* — 2si fc> * c o s < C ) 

It's kind of a modified Pythagorean theorem with a "fudge- 
factor". Note that if C=90 degrees, we da. get the Pythagorean 
theorem.( cos(90*)=0 ) 

Proof: 

A 



c 3 = k 3 + a" 3 and sin(C)=k/b and cos(C)=a'/b 
c 3 = b 3 sin 3 (C) + (a-a') 3 
c 3 = b 3 sin 3 (C) + (a-bcos(C)) 3 
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c 2 = b 2 sin 2 (C) + a 2 -2abcos(C) + b 2 cos 2 (C) 
c 2 = a 2 + b 2 (sin 2 (C)+cos 2 (C)) -2abcos(C) 
c 2 = a 2 + b 2 - 2ab*cos(C) 

In using this theorem, the actual letters are not important, but 
the relative position of the sides and the angle is the important 
thing. In that identity, the angle C and the side c are "opposite 
each other, and the a and b are the other two sides. Thus 
a 2 = b 2 + c 2 -2bc*cos(A) would be an equally valid 

relationship. 

Well, I'll end this "review"(?). You must remember that this is 
traditionally at least a semester high-school course. I've only 
tried to hit the high points as it relates to this book. To 
become competent in.the use of trig, practice with a trig textbook 
is advised. 

V 

GRAPHS OF THE TRIG FUNCTIONS 

"Crank up" your graphing program from the previous chapter and 

insert the SIN and COS functions and see what their graphs look 
_ • _ 
like. .1 might suggest a window of [-12,12]x.[-2,2i to see a good 

representative sample of the functions. Your graphs should look 



1. Find the following:(use the calculator or computer along with 
the definitions given earlier) 

(a) tan(67 degrees) (b) sec(35 degrees) 

(c) esc(1.44) (d) cot(.44) (e) sec(-2.4) 

2. Find the length of side c on the figure to the below. 

^^ 30 * 
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3. Find a formula for these expressions:(hint: 2x=x+x ) 

(a) sln(2x) (b) cos(2x) (c) tan(2x) 

4. Show that cos 2 (x) - sin a (x) = 2cos a (x) - 1. 

5. If you know the sin(A)=1/2 and you know angle A is between 0 
and 90*, what is the value of cos(A)? What if the angle is 
between 90 and 180 degrees? 

6. If you know the sin(A) = 1/2 and the cos(B) = 2/3, what is the 
sin(A+B) and cos(A+B)?(hint: Figure out what cos(A) and sin(B) 
must be and use formulas for sine and cosine of sum of angles as 
given on page 8.) 
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chapter s 1 >c 

PARAMETRIC AND E? OLAR EQUATIONS 



The following topics will give you a more in-depth look at 
different ways that curves can be generated from functions. We 
will first look closer at the general concept of parametric 
equations. 

PARAMETRIC EQUATIONS 

We've already met these when we looked at an alternate form for 
an equation of a line. We'll now generalize the concept and look 
at some different parametric equations. 

DEFINITION: Let I=[a,b]. Let f and g be Real-valued functions on 

I. A curve is a function F:I—>R 2 defined by F((t)=(f(t),g(t)), 
tel. x=f(t) and y=g(t) are called the. parametric equations 
of the curve. t is called the parameter. 

Earlier, parametric equations of a line between two points 
(xi ,yx ) and (x s ,y=i ) were developed. Namely: 

x=xi +(x 2 -xi )t and y=yx +(y= -y= )t, te [0,1]. 

In trigonometry the following picture is 
a common sight:(See to the side) 
x=cos(t), y=sin(t) t e [0,2m] can 
be considered to be parametric equations 
of the unit circle( centered at (0,0) ). 

THEOREM: x=rcos(t)+h and y=rsin(t)+k, t e [0,2m]. can be 

considered to be parametric equations pf a circle of radius r and 
center (h,k). 

Proof: x=rcos(t)+h and y=rsin(t)+k 

x-h=rcos(t) and y-k=rsin(t) 

(x-h) 2 + (y-k) 2 = r 2 cos 2 (t) + r 2 sin 2 (t) 

^ (x-h) 2 + (y-k) 2 = r 2 (cos 2 (t) + sin 2 (t)) 

( (x-h) 2 + (y-k) 2 = r 2 ## 

Often times different curves (different parametric equations) 

result in the same figure, ( or a part of a standard figure ) but 
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the figure is traced out in a different manner. 

Example: x=sin(t) and y=cos(t) / t e [0,2 k] will still result in a 

circle ( x 2 + y 2 =1 ), but the “starting point" is (0,1) (at the 



x=cos(t), y=sin(t) 



Example: x=l/j|r1+t 2 , y=t/S' l+t 3 is the right-hand half of a 

circle. 

Proof: x*+y*= l/(l/t a ) + t^/d+t 2 ) = (l+t^J/d+t*) = 1. 

Since x=l/4^ l+t 52 >0 and y=t/>/ l+t 2 is greater than 0 if t>0 and 
less than 0 if t<0 we only get points in the 1st and 4th 
quadrants. ## 

Convince yourself that no matter how big t gets, the half circle 
never quite is traced out. 

Can you ‘think of other parametric equations that give a curve 
that lies on a circle? 
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(x-h) 2 (y-k) 2 

-- +-= i 

a 2 b 2 . 

A standard parabola with its vertex at (0,0) where c is the 
distance to the focal point is given either by x 2 =4cy or 
y 2 =4cx. 


x 2 =4cy y* =4cx 

If these same parabolas are moved so that their vertex is at 
(h,k) then the equations are: 


(x-h) 2 =4c(y-k) (y-k) 2 =4c(x-h) 

When either of these equations are multiplied out, they are 


easily recognizable 

they are parabolas by the 

fact that 

one 

variable 

is "squared 

" whereas the other one 

isn' 

t. 

- • 


Example: 

y 2 +3x+6 

=0 is a parabola(It's 

not 

in 

standard 

form 

however!) 

There are algebraic tricks to 

put 

such equtions 

in 

standard 

form, but 

I won't get into that 

here. 

( 

By the way 

the 

standard 

form for 

that above equation is 

y a 

= 

4(-3/4)(x+2)— 

vertex is 

at (-2,0) 

and the focal point is 

-3/4 

a 

unit to the 

left 


of the vertex.) 

Even though we can't spend the time going into all the aspects 
of graphing and forming equations of these curves(this is not an 
analytic geometry course), you should take the time to make up a 
few equations, pick a few values for x( or for y) and solve for 
the other variable and plot the points. 
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A usual parametric equation of an ellipse is 

x=acos(t) and y=bsin(t), t [ 0 , 211 ], ( a,b>0 ) 

To show this gives the standard equation, look at the following: 

x=acos(t) and y=bsin(t) 
x/a=cos(t) and y/b=sin(t) 

(x/a) 3 =cos 3 (t) - and (y/b) z =sin z (t) 

(x/a) 3 + (y/b) 2 = cos 2 (t) + sin 2 (t) 

(x/a) 2 + (y/b) z = 1 

This is an equation of a standard ellipse centered at the 
origin. Since cos(t) and sin(t) go between -1 and 1 as t goes 
once "around the circle", the x and y values go between -a and a 
and -b and b respectively. 

As far as I know there is no "usual" parametric equation of a 
parabola. There are several that give parabolas: (a) x=t, y=t 2 
or (b) x=t 52 , y=t or (c) x=sin(t), y=cos(2t) or .... 

To see this last one, recall the trig identity: 

cos(x+y)=cos(x)cos(y)-sin(x)sin(y) 

If you replace both x and y with t, you get 

cos(2t)=cos :2 (t)-sin a (t). Then using the fact that 
sin 2 (t)+cos 3 (t.)=l which is equivalent to cos 2 (t)=l-sin 2 (t), we 
get cos(2t)=l-sin S2 (t)-sin :2 (t) = l-2sin = (t) 

So, y = cos(2t) = l-2sin 2 (t) = l-2x 2 . 

y=l-2x 2 is an equation of a parabola( not 
in standard form however). We are 
describing a curve that is only a part of 
a parabola since -1< sin(t)=x <1 and 
-1< cos(2t)=y <1. See the picture to the 
side of the page for a picture of this 
curve. 

One reason a person doing graphics might wish to be aware of 
different ways to draw a particular curve, is that often the 
visual effect of seeing a curve drawn out is just as important as 
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the curve itself. The process of a curve being drawn is one of 
the most elementary forms of animation. 

pr oblems 

Describe the figure these curves lie on. See if you can change 
these parametric equations to a "normal" x-y equation. Sketch the 
curves, describe the motion of a "bug" traveling along the 
curve(think of t representing time). 

1. x=2t+l, y=3t-4, te[l,5] 

2. x=cos(t 2 ), y=sin(t 2 ), t>0 

3. x=cos(t), y=cos(2t) 

4. x=sin(t), y=sin(2t). Just give its equation in terms of x and 
y. It doesn't have a common name that I'm aware of. Sketch the 
curve for t e [0,2-rc]. 

5. x=|t| , y=t ffi 

6. x=|cos(t)j, y=jsin(t)j 

7. x=tcos(t), y=tsin(t). Don't try to put this in to an x-y 
equation. Just sketch the curve for t>0. 

8. Modify your graphing program to graph parametric equations. 
It shouldn r t take much modification. Your graphing loop will 
depend on T, rather than X. You'll have to ask for input of the 
beginning and ending values of t. 

POLAR EQUATIONS 

In the "real world" there are many ways to describe the position 
of a point in the plane. The most common way is to use the 
Cartesian Coordinate System( Rectangular Coordinate System ). 
This is what has been .used during our discussion of graphing so 
far in this book. Probably the second most common method is the 
Polar Coordinate System. 

Look at the picture to the side of the 
page. The point P is referenced with 
respect to a fixed point 0 and a "half" 


P(r,4) 
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line from that point. It is usual to 
think of 0 to be the same as the origin 
in the rectangular coordinate system and 
the half line to be the positive x-axis. 
This way it is easy to compare 
coordinates of the same point using the 
two systems. Usually the r coordinate 
is a positive distance one moves after 
the angle has been determined. A 
positive angle gives a counter-clockwise 
rotation and a negative angle gives a 
clockwise motion. 

If r is allowed to be negative, the 
convention is to consider the point to be 
r units on the ray that has been rotated 
9 + k . This has the effect that for a 
negative value of r, go "backwards” from 
what you do if r was positive. One 
becomes immediately aware that there are 
MANY( infinite, in fact ) ways to 
represent a point in polar coordinates. 
Consider a few, as you look at the 
diagram to the side. The coordinates 
( 3,te/4 ) or ( 3,9tc/4) or (-3, 5tc/4 ) or 
(-3,-3tc/4) or..... all represent the same 
point. 


P(4, 






This non-unique representation of 
points in the plane is a drawback of this 
particular coordinate system. However, 
the advantages when used in the right 
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situation completely outweighs that one 
disadvantage. 

In the polar coordinate system, certain figures have extremely 
simple equations. Take a circle of radius 2, centered at 0, for 
example. The equation is simply 
r=2. (ie. { (r,8):r=2 } 

A straight line through the origin, oriented 30 degrees from the 
horizontal, has an equation 8 = u/6. (ie. { (r,0):0 =tc/6 } 

It is often the case that certain figures might be easy to 
represent in polar coordinates, but hard to represent in 
rectangular coordinates, and vice-versa. 

You will find that most equations in the polar coordinate system 
are written in the form: 

r = f (0) 

such as: r=sin(28), r=sin(8)+cos(8), r=l-2cos(0), r=0 2 cos(8) 

Often the trig. functions are incorporated as you can see from 
the -examples. 

- Some very fascinating designs can be 
drawn using polar equations. For 
example, r=sin(28) is a 4-leaf rose. 

( 0 e 10,2tc] ) 

Fortunately, for us computer users 
whose screen is basically set up on the 
rectangular coordinate system, there is 
an easy way to convert from polar 
coordinates to rectangular coordinates. 

Basic trigonometry gives us that x/r=cos(0) and y/r=sin(0) 
or x=rcos(8) and y=rsin(0). 

Now, even though in trig we always consider r to positive, if. r 
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is negative the equations x=rcos(0), y=rsin(0) will still give us 
the right x and y coordinates. This is not too easy to verify in 
general, but examine this specific example: 

(4,tc/ 6) or (-4,7 tc/ 6 ) 
x=4cos( tc/ 6)=3.464 y=4sin( tc/6)=2 

or 

x=-4cos(7 tc/ 6)=3.464 y=-4sin(7 tc/6)=2 

It turns out that if we consider r to be negative , we have to 

consider the angle 0 +tc. The values of the sine and cosine 

functions of the different angle will compensate for the fact we 

are using a negative value for r. Note the following: 

-rcos(0+TC)= -r(cos(O)cos(tt) - sin(0)sin(it)) = 

-r(cos(0)(-1) - sin(0)(0) ) = 

-r(-cos(0) - 
rcos(0) 

Similarly -rsin(0+ir)=rsin(0) . I'll leave that as an exercise. 

So, if you've got the program "up and going" for graphing in 
rectangular coordinates, it can easily be converted to graph polar 
equations of the form r=f(0).' Do your FOR-NEXT loop on the 

values of 0 , use your function to calculate a value for r, 

convert to rectangular coordinates by the equations above, and 
plot them as in your regular graphing program. 

A simple outline of such a program might go as follows: 

1. Enter viewport( or default to screen dimensions) 

2. Enter window 

3. Draw axis(optional) 

4. Input beginning and ending values of 0. (Of course you'll have 
to use a different symbol in the program!) 

5. In a FOR-NEXT loop, compute values of r, convert to 
rectangular coordinates, check to see if point is in window, 
convert to screen coordinates, PLOT the point. 
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As one final note, when you consider the equations that convert 
polar coordinates to rectangular coordinates, we can in essence 

think of polar equations being parametric equations with 9 as the 
parameter!! 

pr ob 1 ems 

1. Change the following polar equations to an equation in 
rectangular coordinates. Besides the two conversion formulas 
mentioned in the previous material, you might want to use 

x* + y 2 = r a and/or tan(9)=y/x. 

(a) r=2sin(9) (HINT: mult. both sides by r) 

(b) r=2sin(9)+4cos(9) 

(c) tan(9)=4 (d) r a =sec(29) 

(e) r=2sec(9) (f) r=l/(1-cos(9)) 

2. Sketch the graphs of some of these equations in problem 1 or 
some of the other equations in this chapter. You may want to use 
a graphing program to help you sketch the results. 

3. Convert these polar coord, points to rectangular coordinates, 

(a) (3,30* ) (b) (-4,57* ) (c). (2 -,tc/ 5) (d) (-3.2,-.47) ' 

4. Convert these rectangular coord, to polar coordinates. Give 
two equivalent polar coordinates for each. 

(a) (4,7) (b) (-4,8) (c) (-4,-7) 

5. Finish the proof that -rsin(9+Tt)=rsin(9) 
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CIRCLES 

I shall finish out this chapter with a discussion on graphing 
circles. Circles are common objects to be drawn on a computer 
screen. They are of some interest to us because there are 2 or 3 
distinct ways to draw circles. Also, they are of interest because 
typically circles are slow to draw out because of the calculations 
that must be performed in most drawing routines. It is then of 
interest to see if one technique is faster than another. 

THE OLD STANDBY 

(x-h) 3 + (y-k) 3 = r 3 is the standard 

equation of a circle of radius r centered 
at <h,k). Or, solving for y, we get 
y=±\/ r 3 - (x-h) 3 + k 

Since the computer can't handle this "± jazz", we'd have to look 
at this as two functions: 

y=vr 3 -(x-h) 3 +k y= —V r 3 -(x-h) 3 +k 

A sample routine might look like this: 

( I'll not use the window routines in these short examples, hence 

the circles will graph themselves out "upside down" . You could 

consider the window to equal the viewport which is 

[0,639]x[0,199]. ) 

LET r=30 ’ 

LET h=100 
LET k=60 

FOR x=h-r TO h+r 
LET y=SQR(r A 2-(x-h) A 2) 

PLOT x,y+k 
PLOT x,-y+k 
NEXT x 

Just in this little example, we can experiment with the "speed" 
concept. Try replacing the first line in the loop with 
y=SQR(r*r-(x-h)*(x-h)) . I think you will see that simple 
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multiplication is performed much faster than exponentation. 

A more common way to view the circle at (h,k) is to think of it 
as a circle at (0,0) that's been moved to (h,k). In other words 
use x 2 + y a =r a . 

(or y=±i/r 2 -x = )7 then plot (x+h,y+k). 

LET r=3Q 
LET h=100 
LET k=60 
FOR x=-r TO r 

LET y=SBR(r*r-x*x) ! X*X IS FASTER THAN X A 2 
PLOT x+h,y+k ! TRANSLATES THE POINT BY AN AHOUNT H & K 
PLOT x+h,-y+k ! DITTO 
NEXT x 

One could hasten the graphing by making use of more symmetry. 

Note the picture to the right: 

LET r=30 
LET h=10Q 
LET k=60 
FOR x=-r TO r 
LET y=SBR(r*r-x*x) 

PLOT x+h,y+k 
PLOT x+h.-y+k 
PLOT -x+h,y+k 
PLOT -x+h,-y+k 
NEXT x 

Adding h or k moved our points to where we wanted them. This 
concept of translating a picture to someplace else will be a major 
topic in a later chapter. 

Can you write a program that draws a circle that takes even more 
advantage of the symmetry of the circle? 


TRIG APPROACHES 

From our previous discussion of the 
basis of trigonometry, we are familiar 
with the diagram to the right. Thus for 


a circle at the origin (0,0) the 
parametric equations 


x=rcos(t) and y=rsin(t) , 05ti2n 
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will do the job. To move the circle to a new origin (h,k) simply 

add h to x and k to y. 

LET r=30 
LET h=100 
LET y=60 

FOR t=0 TO 6.3 STEP 3.14/60 ! 3 DEGREE INCREMENTS 
LET x=r*C0S(t) 

LET y=r*SIN(t) 

PLOT x+h,y+k 
NEXT t 

(or replace the loop with:) 


FOR t=0 TO 1.57 STEP 3.14/60 ! 0 TO 90 DEG STEP 3 DEG 
LET x=r*C0S(t) 

LET t=r*SIN(t) 

PLOT x+h,y+k 
PLOT -x+h,y+k 
PLOT x+h,-y+k 
PLOT -x+hj-y+k 
NEXT t 

(Or replace the loop with:) 

PLOT r*C0S(0)+h,r*SIN(0)+k ! PLOT FIRST POINT ON CIRCLE 

FOR t=0 TO 3.14 STEP 3.14/60 ! 0 TO 360 DEG STEP 3 DEG 
LET x=r*C0S(t) 

LET y=r*SIN(t) 

DRAW TO x+h,y+k ! DRAW LINE TO NEXT POINT 
NEXT t 

Trig routines can be fairly slow-because it is very time 
consuming to calculate values of trig functions. 


A FASTER WAY 

The two following formulas from trig will guide us toward more 
speed: 

COS(A+B) =COS(A)COS(B)-SIN(A)SIN(B) 

SIN(A+B)= SIN(A)COS(B)+COS(A)SIN(B) 

Consider the diagram to the right: 

Suppose we're trying to go from point 
(Xr»,yr. ) to the next point (x n +i,yn+i ), 
incrementing the angle by a fixed amount 
dt. Using trig definitions for 
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coordinates of a point on the circle we 
get: 

Xn = rcos(t) yr, = rsin(t) 

Xn+i = rcos(t+dt) y n +i = rsin(t+dt) 

Let's now play with x„ + i and y n+ i by using our trig formulas. 
Xn+i = rcos(t+dt) = r[ cos(t)cos(dt)-sin(t)sin(dt) ] 

= rcos(t)cos(dt)-rsin(t)sin(dt) 

= Xn cos(dt) - y n sin(dt) 


Similarly 

yn-i = rsin(t+dt) = you fill in 

= the blanks 

= yn cos(dt) + Xn sin(dt) 

Since cos(dt) and sin(dt) would simply be constants, we easily 

get new points on our circle by using the previous point plotted 

and multiplying by constants. 

LET r=30 
LET h=10Q 
LET k=60 

LET dt=3.14/60 ! A 3 degree increaent 
LET s=SIN(dt) ! A tise-saver 
LET c=C0S(dt) 

LET x=r ! beginning value of x 

LET y=0 ! beginning value of y 

FOR t=0 TO 6.3 STEP dt 
PLOT x+h,y+k 

LET xl=x*c-y*s !Now coapute the next point 
LET yl=y*c+x*s 
LET x=xl 
LET y=yl 
NEXT t 

The above routine draws out a nice looking circle clockwise 
(because . our screen is upside down) in a very efficient manner. 
For further speed( sacrificing some elegance in animation) use the 
symmetry as outlined in the preceeding programs. 

I'll finish by including(without comment) a very fast circle 
drawing routine. Note that there is no trig, square roots or even 
multiplcation. 
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LET r=30 
LET h=100 
LET k=60 
LET phi=0 
LET x=r 
LET y=0 
REPEAT 

PLOT x+h,y+k — 

PLOT -x+h,y+k 
PLOT x+h,-y+k 
PLOT -x+h,-y+k 
PLOT y+h,x+k 
PLOT -y+h,x+k 
PLOT y+h,-x+k 
PLOT -y+h,-x+k 
LET phiy=phi+y+y+i 
LET phixy=phiy-x-x+l 
LET phi=phiy 
INC y 

IF ABS(phixyXABSCphiy) THEN 
phi=phixy 
DEC x 
ENDIF 
UNTIL x<y 

Why don't you see if you can put all of these programs together 
as one program to serve as a Circle demo program? I think you'd 
find it interesting to see how the circles draw out, one right 
after the other. Of course adjust the values of H and K so the 
circles don't all draw on top of each other. 



ST Math & Graphics 


chapter seven 


CL I F>E> I NG 

When condsidering a window (and its corresponding viewport) in 
doing graphics, it becomes important to determine what's in the 
window and what li-es outside the window. 

In displaying a complex shaped object, such as the one pictured 
to the right, several types o£ cases arise: 

(1) Line segments whose endpoints are 
entirely contained in the window (line 


FG). 

(2) Line segments whose endpoints are 
not in the window, but part of the line 
intersects the window (line AB). 

(3) Line segments whose endpoints are 
not in the window and no part of the line 
is in the window (line DE). 


B 



(4) Line segments where one endpoint is 
in the window and other endpoint is not 
in the window (lines BC,CD,EF,AG). 

Displaying only the portions of the object that intersect the 
window is called CLIPPING. It becomes the task in this chapter to 
develop an alogrithm and program that will perform this clipping 
process. 


FIRST ALGORITHM 

If one is graphing the lines, point by point, by using one of 
the earlier routines, the job is simple: compare the (x,y) 
coordinate of the point to be plotted with the minimum and maximum 
x and y values of the window. If the point lies in the window, 
plot it, otherwise don’t plot it. This is often an unsatisfactory 
solution when working with a high-level language due to the 
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slowness of line plotting by plotting it point by point. However, 
in machine language this algorithm is extremely practical! 

SECOND ALGORITHM 

The algorithm we will develop assumes 
you will use the "built-in" drawing 
capability of the language. 

(ie. DRAW xl,yl TO x2,y2 ). 

The algorithm, simply stated, is this: 
if one or both endpoints lie outside the 
window, determine the intersection points 
with the window. Then draw the line 
between the two points in (or on) the 
window. 




The problem with this algorithm is the 
difficulty in writing a routine that 
covers all the possible positions of the 
line segments. There are so many 
possible orientations of the two points, 
and seemingly small changes result in 
different types of intersections with the 
window. 



Sometimes you have no intersections; sometimes 1 intersection, 
either with a horizontal side or a vertical side; sometimes 2 
intersections, both with vertical sides, or perhaps both with 
horizontal sides, or one with a vertical side and one with a 
horizontal side. Trying to write a routine that takes into 
account all of these possibilites, and yet.is fairly quick and 
efficient becomes an interesting task. 

The heart of the routine that will be 
shown you is the use of the parametric 
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equations of a line: 

X=X * +(X 2 -X 1 )T and Y=Yx+(Y 2 -Yi)T 
where when T=0, we get (Xi,Yi) and when 
T=l, we get (Xs^Ya), and for any value of 
T between 0 and 1 we get points on the 
line segment between these two points. 

We will also be interested in these equations solved for T: 

T=(X-Xi)/DX and T=(Y-Yi)/DY 
where DX=(Xa-Xt) and DY=(Ya-Y ±) 

The program that follows (at end of chapter) was written by a 
student of mine the first time I taught a class in Math and 
Graphics. I had presented my routine, and he saw a better way to 
do the job. I invite you to build a better "mousetrap"! 

I won't explain every part of the program as the ideas are 
duplicated for the different situations. Examine the listing that 
follows as the main parts are commented upon. 

LINES beginning with IF (yl<ys AND y2<ys) OR. (........ . 

This -checks to see if both points are 
above, below, or to the sides of the 
viewport (The boundaries of the viewport 
are given by the variables xs,xl,ys,yl). 

If they are as shown in the figure to the 
right, then the subroutine RETURNS to 
main program witthout drawing a line. 

LINES beginning with LET FLAG=0 . 

The routine assumes the vector is going from left to right, with 
(xl,yl) being the left-end point. The points are swapped if 
necessary. A variable called flag is set to 1 if a swap is made. 
The flag is checked before leaving this procedure(See last few 
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lines of this subroutine) and the points are "unswapped" if they 
had been swapped originally. 

LINES beginning with Let dx=x2-xl ...... 

The values of x2-xl and y2-yl are used repeatedly. Call these 
dx and dy respectively. The program then performs one of three 
actions if the line is vertical (dx=0), depending on whether dy is 
negulive, zero, or positive. 


LINES beginning with ELSE ! dx is not 0 ....... 

The program performs 1 of 3 actions depending if the line is 
uphill (dy<0—backwards from normal since the y-axis on the screen 
is upside-down), horizontal (dy=0), or downhill (dy>0). 

Let's consider the THIRD "If 
statement"(If dy>0 THEN ). The lines have 
an orientation as depicted to the right. 

Our explanation of what happens is 
typical of the process in the other "If 



statements". 

Be sure to examine the figure to .the 
right as this explanation proceeds. The 
figure shows the lines forming the 

viewport and the line segment itself 
being extended to show "intersections". 
Remember the line is being drawn left to 
right. A line that crosses the viewport 
will- intersect the left side or top 
before it would hit either a right side 
or bottom (convince yourself this is not 
true for a line that misses the 

viewport). 
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To find the values of t for the left or top substitute xs or ys 
in for x or y in the parametric equations that were solved for 
t.(See top of page 3) Call one of the values t9 and the other tb. 
tb would be the larger of the two values( tb=MAX(tb,t9) ). If tb 
is less than 0 then our starting point is within the boundaries of 
the viewport, so let tb=0(Recall that t=0 gives (xl,yl)—our 
starting point). Now we check the intersection with the right and 
bolLom by letting x be xl and y be yl. We'll call these values t9 
and te. We want te to be the smallest of these two values 
(te=MIN(te,t9) ). If te>l then the point was in the window, so we 
want te=l. Note in the picture to the right that "normally" tb is 
less than te if the line crosses the viewport. If tb is greater 
than te then the line misses the viewport. If te and tb are ok, 
then the x and y values associated with those values of t are 
calculated at the spot labeled draw_line at the end of the 
program. Then a line is drawn from (xb,yb) to (xe,ye). 

It is hoped that you will examine the rest of the listing 
to thoroughly understand the logic. . Don't be discouraged if you 
have trouble though. I found it hard to understand other clipping 
routines myself and finally wrote my own(which I understand!). I 
invite you to write your own routine also. 

If you write your own routine or type in mine, be sure to save it 
as an ASCII file (SAVE,A option) so you can MERGE it with your 
other programs. 

DEMO #l:(At end of chapter.) 

In this little program, a series of diagonal lines are created 
in a FOR-NEXT loop. The program draws a viewport, then through 
the use of the clipping subroutine, draws only the portion of the 
lines that go through the viewport. If you type in the demo, be 
sure Lo MERGE in the clip_line procedure. 



ST Math & Graphics 


DEMO #2:(At end of chapter.) 

This program draws a picture that is described in DATA 
statements. The figure is described in a window of 
[0,100]x[0,100]. Thus all the coordinates are positive. I take 
advantage of this and use negative numbers as flags in my Data 
lines.You can display this window or any part of that window on 
any viewport you wish. The program will draw the part of the 
figure that belongs in the viewport. In this program, the lines 
are drawn sequentially, hence the second point for the first line, 
becomes the first point of the next line to be drawn. 

Be sure to try to find my initials on the door by picking the 
right window! 

Be sure to look at the DATA lines to see how certain numbers in 
the data are used as flags to indicate that something else is 
happening besides drawing points. That's why those -l's and -9's 
are in the data list. The -l's indicate the pen is to be lifted 
and start a new sequence of lines. The -9's indicate the end of 
the data list.. 

pr oga mm ± ra <g project 

Use the Demo #2 program as a model and draw your own picture. 

Include some real small detail that won't show up unless the user 
picks a window that will magnify that portion of your picture. 
You might try a couple of levels of "smallness". 

You might try to go beyound this example and try to add some 
color to your drawings. Try to set flags in your DATA that 
indicates you wish to change color (similar to the way you 
indicated to the program you wished to "lift your pencil" to start 
drawing in a new place. 
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PROCEDURE dip.line 

' clipping subroutine--by Ken Hall—Highline CC--'84 
1 X5,xl,ys,yl--viewport diaensions 
' xl,yl,c2,y2—endpoints of line 
LOCAL t9 f tb,te,xb f yb,xe f ye f flag 

1 tb,yb,xb are values associated with the beginning point 
’ te,ye,xe are values associated with the ending point 

f 

IF <yl<ys AND y2<ys) OR (yl>yl AND y2>yl) OR (xl<xs AND x2<xs) OR (xl>xl AND x2>xl) THEN 
GOTO exit sub 
ENDIF 
LET flag=0 
IF x2<xl THEN 
SNAP xl,x2 
SWAP yl,y2 
LET flag=l 
ENDIF 

LET dx=x2-xl 
LET dy=y2-yl 
IF dx=0 THEN 
IF dy<0 THEN 
tb=NAX((yl-yl)/dy,0) 
te=HIN((ys-yi)/dy,1) 

GOTO draw_line 
ENDIF 

IF dy=0 THEN 
PLOT xl,yl 
GOTO exit_sub 
ENDIF 

IF dy>0 THEN 
tb=HAX((ys-yl)/dy,0) 
te=HIN((yl-yi)/dy,l) 

GOTO draw_line 
ENDIF " 

ELSE ! Dx is not 0 
IF dy<Q THEN 
t9=(xs-xl)/dx 
tb=(yl-yl)/dy 
tb=HAXCtS,tb) 
tb=NAX(tb,0) 
t9=(xl-xl)/dx 
te=(ys-yl)/dy 
te=NIN(te,t9) 
te=MIN(te,l) 

IF te<tb THEN 
GOTO exit sub 
ENDIF 

GOTO draw line 
ENDIF 

IF dy=0 THEN 
tb=NAX((xs-xl)/dx,0) 
te=NIN((xl-xl)/dx,l) 

GOTO drau_line 
ENDIF 

IF dy>0 THEN 
t9=(xs-xl)/dx 
tb=(ys-yl)/dy 
tb=HAX(tb,t9) 
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tb=HAX(tb,0) 

t9=(xl-xl)/dx 

te=(yl-yl)/dy 

te=HIN(te,t9) 

te=HIN(te,l) 

IF te<tb THEN 
GOTO exit sub 
ENDIF 

GOTO draw line 
ENDIF 
ENDIF 

drawjine: 
xb=xl+tb*dx 
yb=yl+tb*dy 
xe=xl+te*dx 
ye=yl+te*dy 
DRAH xb,yb TO xe,ye 
’ *«•»«*»»****<**** 
exit_sub: 

IF flag=i THEN 
SWAP xl,x2 
SWAP yl,y2 
ENDIF 
RETURN 

JHHHmHHHlHHHHHHiHHHHmHTOHHHWHHH 
' Dbibo il of Clip_line subroutine 
INPUT "enter viewport dis.(xs,xl,ys,yl) ”,xs f xl,ys,yl 
BOX xs,ys,xl,yl 
FOR z=-140 TO 200 STEP 3 
xl=0 

yl=140+z 

x2=500 

y2=r 

clip line 
NEXT I - 

’ HERE IS WHERE TO HERGE IN YOUR CLIPLINE SUBROUTINE 
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' Deao #2 of using Clipjine routine 
' Deso of using clipping routine 
’ ft house is drawn using data statements 

i 

' try window of 60,70,18,25 to see initials on door 

INPUT "Enter viewport dimensionsfie._5,400,4,150) ",xs,xl,ys,yl 

INPUT "enter windowfie. 0,100,0,100 ) ",a,b,c,d 

BOX xs,ys,xl,yl 

LET wi=b-a 

LET w2=d-c 

LET vl=xl-xs 

LET v2=ys-yl 

RESTORE houedata 

DO 

READ x,y ! pen is up, new line 

EXIT IF x=-3 

window_function 

LET xl=x 

LET yl=y 

DO 

READ x,y 

EXIT IF x=-l OR x=-3 
window_function 
LET x2=x 
LET y2=y 
clipjine 
LET xl=x2 
LET y1=y2 
LOOP 

EXIT IF x=-9 
LOOP 

hoaedata: ‘ 

DATA 20,50,20,10,80,10,90,50,20,50,30,70,70, 50,20,50,-1,-1 / . - 

DATA 30,30,30,20,50,20,50,30,30,30,-1,-1 , 

DATA 60,10,60,25,70,25,70,10,-1,-1 

DATA 45,70,45,85,55,85,55,70,-1,-1 

DATA 45,74,55,74,-1,-1,45,78,55,78,-1,-1,45,82,55,82,-1,-1 

DATA 47,85,47,82,-1,-1,53,85,53,82,-1,-1 

DATA 46,82,46,78,-1,-1,50,82,50,78,-1,-1,54,82,54,78,-1,-1 

DATA 49,78,49,74,-1,-1,51,78,51,74,-1,-1,54,78,54,74,-1,-1 

DATA 47,74,47,70,-1,-1,53,74,53,70,-1,-1 

DATA 62,20,62,22,68,22,68,20,62,20,-1,-1 

DATA 63,21,63,20.3,64.5,20.3,64.5,21,63,21,63,21.7,64,2,21.7,64.2,21,-1,-1 
DATA 65.5,20.3,65.5,21.7,-1,-1,65,5,21,67,21,-1,-1,67,21.7,67,20,3 
DATA -3,-9 

PROCEDURE vindowjunction 
x=(x-a)/wl*vl+xs 
y=(y-c)/w2*v2+yl 
RETURN 

' HERE IS WHERE TO MERGE IN YOUR CLIPLINE SUBROUTINE 

********************************************************** 
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chapter eight 
MA.TR. I CE S 

Before we can go much further in the traditional topics in 
graphics, I need to acquaint you with a new(?) concept. This 
chapter will essentially be "pure" math—later chapters will 
connect the topic to graphics. 

MATRICES 

A matrix is a rectangular array of numbers. An mxn matrix is a 
matrix with m rows and n columns. 

Examples: 

|2 3 4| 

is a 2x3 matrix 


A= 


2 3 

3 7 


B= 


C= 


3 4 

4 9 

2 4 

2 1 
4 7 

2 8 

3 3 


is a 3x3 matrix 


is a 4x2 matrix 


D= [3 4 7 9 8] is a 1x5 matrix 

A square matrix is a matrix with the same number of rows and 
columns. Matrix B above is a square matrix. 


A matrix of this type 


10 0 0 
0 10 0 
0 0 1 0 
0 0 0 1 


is called an identity matrix 


8-1 



ST Math & Graphics 

The following notation is often us.ed when writing a matrix, 
using letters to represent the entries in the matrix. 


dll dl2 dl3 ••••••• » din 

dizzx 3 22 3 . 23 . .. .d2n 

331. d-32 333 .. 33n 

A= 


|3m± 3 m2 3 m3 ••••«•« • 3mr»J 

This subscript notation allows one to refer to matrices in 
general or to make definitions of particular matrices. For 
example an identity matrix I can be described as follows: 

I = ta±j] where fa*!. = 1 if i=l,2,..n 

\a±j =0 if i does not equal jj 

MATRICES IN BASIC 

Most versions of BASIC allow you to establish a matrix as a 
means of storing numbers. 

In a program, a line like DIM A(7,5) will reserve memory 
space for a matrix of & rows and 6 columns (. or .7 rows and 5 
columns in some BASICS.) Generally there is a row "zero" and a 
column "zero". The DIM A(7,5) statement establishes the largest 
row and column number, so depending on whether your computer 
allows a 0 row and column, you get either 7 rows, 5 columns or 8 
rows, 6 columns. 

A statement in your program like Let A(4,3)=5 puts the 
number 5 in the spot that is referenced by row 4, column 3. 

A statement in your program like Let T=A(3,2) assigns to the 
variable T the value found in the matrix A at the spot referenced 
by row 3 and column 2, 

Hence you can see that in BASIC when we write A(3,2), this is 
equivalent to refering to a 3 » .(Assuming A=[aij] ) 
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WHAT ARE THEY GOOD FOR? 

Matrices might be set up for any number of reasons. Often it is 
just a good way to organize data. 

Examples: 

- Warehouse 



A 

B 

C 

Chairs 

500 

10 

200 

Sofas 

50 

40 

20 

Tables 

300 

30 

0_ 


A matrix filled with l's and 0's can be used to represent 
alphabetical characters to be drawn on a screen. 


r 

The 1 tells the computer to plot a point, 0 means don't plot a 
point.' 

Sometimes the matrix might represent a series of points or vectors 
that are used to draw a figure. 


nt a series of simultaneous equations to 

'3 4 1 7' 

1 - 12-1 
2 10 4 


20 

5 
2 
1 

6 


30 

4 

5 
7 

89 


A matrix might represe 

be solved. 

3x + 4y + z = 7 
x - y + 2z = -1 
2x + y =4 


A crude 



representation of A 
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in this book we will primarily be using matrices to represent 
transformations of points in doing graphics in the plane and 

3-dimensional space. 

MATRIX OPERATIONS 

Besides representing "things", matrices can be made to behave 
like mathematical objects in their own right. To do this we need 
to define some operations that are commonly performed. 


ADDITION 

Let A=[a=Lj] and B=tbi_j] be m x n matrices. 
Define A+B=[a±j +btj ]. 


Example: 

2 3 4 


-1 2 6 

A- 

3 1 7 

B= 

. 3 2 4_ 


A+B= 


1 5 
6 3 


10 

11 


The matrix 0=[01 is called the additive identity matrix since 
it would be true that A+0=A. 

SCALAR MULTIPLICATION 

Let A=[a±j] and let t be any real number(t is called a scalar). 
Define tA=[ta±j]. 

Example: 1*2 3 4] 


A= 1 5 7 

_3 2 1 


Then 

4 6 8 


-6 -9 -12' 

2A= 

2 10 14 

-3A= 

-3 -15 -21 


6 4 2_ 


-9 -6 -3 


These two operations should remind you of our operations for 
vectors. If you think about it, a vector in the plane is a simple 
1x2 matrix. As you will see later, this way of looking at a 
vector will be important. 



ST Math & Graphics 


MATRIX MULTIPLICATION 

Let A=[a±j] be an m x n matrix and B=[b±j] be an n x p matrix. 

n 

AB is the m x p matrix [E an< b*j ]. 

*<— i 

n 

In other words if C=[c ±J ]=AB, then c ±J = E a±„ b^j 

»<— 1 


As you will see after we go through a couple of examples the 
number in row i, column j is obtained by multiplying row i of A 
with column j of B, adding as we go. 


Example: 2 ll 

A= 3 2 

4 1 

3x2 matrix 

AB is the 3x4 matrix: 


B= 


13 7 4 

_2 1 5 2_ 

2x4 matrix 


2 *1+1*2 

2*3+l*l 

2*7 + l*5 

2 *4+1 *2* 


4 

7 

19 

6 “ 

3*l+2*2 

3*3 + 2*l 

3*7+2*5 

3*4 + 2*2 

= 

7 

11 

31 

16 

4*1+1*2 

4*3 + l*l 

4*7 + l*5 

4 * 4 +1 * 2_ 


6 

13 

33 

!8 


This may seem like a complicated procedure, . but it actually is 
pretty easy to do with a little practice. I find my hands really 
getting into the action, using a finger on the left hand to go 
across a row of the first matrix and a finger on my right hand to 
go down a column in the second matrix. It even helps to "mutter” 
out loud!! Look again at our previous example: 


2 1 


1 3 7. 4' 


4 7 19 6 

3 .2 


2 1 5^ 2 

II 

7 11 @ 16 

4 1 



6 13 33 18 


The number 31 which was in row 2, 
"multiplying row 2 by column 3". 


column 3 was obtained by 
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Example: Take our previous warehouse matrix. 


Suppose each chair has a value of $200. 
Suppose each sofa has a value of $500. 
Suppose each table-has a value of $300. 


Arrange in a 1x3 matrix: 



A 

B 

C 

Chairs 

500 

10 

20ol 

Sofas 

50 

40 

20 

Tables 

300 

30 

0 _ 

lir Sofa 

Table 



[200 500 300] 


500 

10 

200 

50 

40 

20 

300 

30 

0 


Value [200 500 300] 


= [215,000 31,000 50,000] 


represents the value of the inventory in each warehouse. 


Example: The system of equations 2x+3y+4z=10 

5x+2y-z =3 is often 
2x+y+3z =4 


represented by 

problems 

h 2 


2 3 4 


X 


10 * 

5 2-1 


y 

= 

3 

2 1 3. 




_ 4, 


in linear algebra 


1. Let A= 


1 4 


B= 


2 1 
3 5 


C= 


\2 3 -1 
4 5 7 


(a) Find A+B; AB; BC; 5B; BA 

(b) Is matrix multiplication commutative?(Do any of your results 

above confirm your answer?) 

\ - - 

(c) Why isn't CB defined in this example? 



'3-4 2 


2.1 2 

2. Let A= 

5 6 2 

B=. 

0 10 


-1 2 -3 

. 


,5 2 1 


AB=? BA=? A a =AA=? 
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3. The matrix -A is defined to be the same as -1A. Using A as 
defined in problem 2, find -A. Is there a "reasonable way" to 
define subtraction of matrices—do so if you think there is. 



3 7 2 


1 o o" 

II 

< 

2 4 1 

1 = 

0 10 


2 1 2 


0 0 1_ 

(a) Find AI and 

IA. 

In word 

the matrix 

I . 



(b) 

i o o' 



Let K= 

0 0 1 

Find AK a 


0 1 0. 




In words summarize your conclusions about the matrix K and its 
effect on other matrices. 


(c) 


Let N= 


-10 0 
0 10 

0 0 1 


Find AN and NA. 


In words summarize your conclusions about the matrix N and its 
effect on other matrices. 

5. Can you describe any other possible uses for matrices?(You 
might think of games, maps, photographs...) 
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chapter n±ne 

MATRICES AND TRANS FORMATIONS 

OF FOINTS 


A square matrix of the form [A C 
transformation of - points B D 

■ - H 

Note: 


can be thought of as a 
in the plane. 


Cx y] 


A 

B 


C 

D. 


[Ax+By Cx+Dy] 


We have transformed a point (x,y) (thought of as a 1x2 matrix) 
to a new point (Ax+By / Cx+Dy). Note that in this case our matrix 
that did the transformation could be thought of as a function from 
R a —>R a . 

It might be noted we could get-the same result by writing: 


A B 


X 


Ax+By 

C D 


y. 


Cx+Dy 


personally, kind of prefer this second form. However in most 
of the literature about, transformations, the first form seems to 
prevail. Just as an aside, note the first matrix and the second 
matrix are related by the fact that their rows and columns have 
been interchanged. If this happens, often the one matrix is 
called the transpose of the other matrix. 


fA C 


For a transformation 


LB DJ 

it will be left unchanged. 


, the origin is special: 
To see this note that 


[0 0 ] 


A 

B 


C 

D. 


=[0*A+0*B 


0*C+0*D]=[0 0] 


* 


Also transformations of this type "preserve" lines. That is, 
points on a line, get transformed to points on another line. To 
see this, imagine we have a line segment between (xi, y ±) and 
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(x 2 ,y 2 ) . (Xx,yx) gets transformed into a point (x t ',yi') 
and (x a/ y 2 ) becomes a point (x 2 ',y 2 '). For example. 


[xi',yi'] 


[Xx 


Ta 


yx ] 


l B 


c 

D. 


[Axx +Byx Cxi +Dyx ] 


Let’s see if a point on a line between (Xx,yx) and (x 2 ; y 2 ) 
becomes a point between (Xx',yx') and (x 2 ',y a '). 

( Xx+(X S -Xx)t, yx+(y*-yOt ), 0<t<l is a point between (x±,yx) 


and (x=>,y s ) . 


[ Xx+(X 2 -Xx)t 


yi+(y2-yi)t ] 


A 

B 


C 

D 


[ Axx+A(X=»-Xx ) t+Byx+B(y = -yx) t Cxi+C( x 2 -Xi ) t +Dyx+D(y = -y x )t ] = 

[(AXx+Byx)+((Ax 2 +By = )-(Axx+Byx))t (Cxx+Dyi)+((Cxa+Dy = )-(CXx+Dyx))t] 
= tx.»’ + (x*'-xx’)t yx' + (y 2 ’-yx')t ] ,0<tsl 
These are, as you can see, the parametric equations of a point 
on the line segment determined by (Xx',y»'> and (x 2 ',y 2 ' ). such 
a transformation is called a linear transformation. 

Since, lines remain lines, this means other figures pretty much 
retain their shapes.( They may get distorted, however.) This is 
because you generally can approximate any curve with "short" line 
segments. 


TYE»ES OB’ TRANSFORMATIONS 
Scaling 

The size of an object may be controlled by a matrix of the 


type: 


Sx 0 

0 S 




[x y] 


= [Sx x S y y] 


S~ 0 
0 Syj 

Note that the x coordinate was multiplied by Sx and the y 
coordinate by S x . This will give a "stretching" effect if Sx or 
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S y >1 amd a "shrinking" effect if S* or are between 0 and- 1. 





If S«=Sy f then overall uniform scaling takes place. 



(toward) the origin. Remember the origin, if it is part of the 
figure you're scaling, does not move. These facts must be given 
consideration when trying for certain effects. Experiment by 
drawing simple objects and trying different combinations of 
positive numbers for S* and s*. . 

Reflection 

These matrices cause symmetric movement of points. 


1 0 | 


About the x-axis: 


l° -y 
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A matrix such as 
reflection. 



0 

3 


is a combination of scaling and 



-2 O’ 


2 0~ 


-1 o' 


■-1 o' 


2 6 

Note that 


~ = 




or 





. 0 3. 


0 3 


. 0 1, 


. 0 1 


? 3 



application of 


cos9 sin8 


0 > 0 


|-sine cosGJ 

( 0 < 0 gives rotation in clockwise direction.) 
To see that this matrix gives the desired result. 



consider 


From trig we know that x=rcos$ and y=rsin$ 

Also, looking at (x',y'): x'=rcos( $ + 0 ) and y'=rsin( $ + 0 


the 
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Applying the identities of cos(a+b)=cos(a)cos(b)-sin(a)sin(b) 
and sin(a+b)=sin(a)cos(b)+cos(a)sin(b) 


We get x =rcos($ +0 )=r(cos$ cose -sin$ sine ) 
=rcos<i> cose - rsin<i> sine 
=xcose - ysine 


Similarly 


y* =rsin($ +0 )=r(sin$ cose +cos$ sine ) 
=rsin , i> cose + rcos$ sine 
=xsin0 + ycos© 


Thus [x' y']=[x y] 


cos© sinB 
-sine cosB 




Note that rotation matrices always assumes the origin is the 
center of rotation. This may cause a surprise if you forget 
this!* 



.876 


Translation 



Translation is a simple concept: simply add a constant to the x 
and y coordinates: x’ =x+h y* =y+k . However, you can not 
describe this transformation using 2x2 matrices. Rather than 
giving up on the idea of using matrices to describe this 
transformation, someone(?) came up with this ingenious idea: 

Define ( x / ¥r 1) = ( x ,y) • This is called the homogeneous 
coordinates of (x,y). 
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In place of our usual 


A C 
B D 


use 


A C 0 
B D 0 
0 0 1 , 


For translation by h and k use 


10 0 
0 10 
h k 1 


[x y]=[x y 1] 


10 0 
0 10 
h k 1 


= [x+h y+k 1] = tx+h y+k] 


When combining transformations, deciding on which order things 
occur can be important. Look at the figures on the previous page. 

Suppose we take as our original object the square with the lower 
left point at the origin. The top figures show a rotation of 30 
degrees. Now imagine a translation of h=l, k=l (Not shown—sketch 
it in.) The bottom set of figures shows a translation of h=l, k=l 
and then a rotation of 30 degrees. You should convince yourself 
that the result is not the same!! 

The programs that follow utilize these transformation fo'rmulas. 

The first is a simple demo of transforming a simple object. You 
can not do multiple transformations (ie. progressively 
transforming an object). 

The second program utilizes all of the "bells and whistles"— 
window, viewport, clipping, progressive transformations. 

The third program is an "application" program that utilizes the 
ideas to draw "snowflakes". 

The fourth program allows you to experiment with the 

transformation matrix in showing a familiar Northwest figure. 


More about the 1st program: 

This program assumes the viewport is the entire screen. 
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xC and yC are the center coordinates o£ the screen. 

The window is basically t-xl/2,xl/2]x[-yl/2,yl/2] 

The mapping formulas from the window to the viewport becomes .very 

simple: x=x+xc;y=-y+yc 

The program uses "ho clipping. 

The program has the following simple outline: 

(1) Initialize viewport, window, center of screen 

(2) Draw x-y axis 

(3) Get desired transformation 

(4) Read Data, applying transformation 

(5) Change to screen coordinates 

(6) Draw the figure 
problems 

1. Give the matrix that would cause a vertical scaling of 5 and a 
horizontal scaling of 2. 

2. Give the matrix that would cause an overall reduction in the 
size of a figure by a factor of 1/2. 

3. Give the matrix that would cause a rotation of 90 degrees. 
Give a matrix that would cause a rotation of -45 degrees. 

(Load up the last program(Seahawk) from your disk, and test your 
answers in problems 1-3—to see the original seahawk, make the 
matirx be the Scaling matrix with S* and S y both be 1.) 

4. Give a matrix that would cause a translation of a figure 10 
units to the right and 25 units down. 

5. Give a matrix that would give the effect of a translation to 
the right 20 units, then a reflection about the y-axis.[Hint: 
consider multiplying matrices together.] Give a matrix that would 
reflect about the y-axis, then translate to the right by 20 
units. Are these the same matrices? 
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’ siaple deao of transforaations 
’ Siaptran 

' Medius Resolution 
LET xl=639 
LET yl=153 

LET xc=INT(xl/2) CENTER OF SCREEN 
LET yc=INT(yl/2) 

' DRAH IN AXIS 

DRAH 0,yc TO xl,yc ! Y AXIS 

PRINT AT(77,12);INT(2*yc) 

DRAW xc,0 TO xc,yl ! X AXIS 
PRINT AT(42,l)jINT(yc) 

read_and_plot('H“) ! plot the original figure (H=K=0) 

DO 

' TRANSLATION MENU 
REPEAT 

PRINT AT(1,21)j"Choose Move Scale Rotate Suit (M,S,R 
or 3)>>“| 

LET choice*=UPPER*(INPUT$(1)) 

UNTIL INSTR("MSRS",choice*)>0 
PRINT choice* 

EXIT IF choice*=“B“ 

IF choice$=“M" THEN 

INPUT “Enter your translation-- H,K (return) »’,h,k 
ENDIF . 

IF choice*=“S" THEN 

INPUT “Enter the scale factors-- Sx,Sy (return) »",sx l sy 
ENDIF 

IF choice$="R" THEN 

INPUT “Enter your rot at ion (degrees)--Theta (return) »“,theta 
LET theta=theta*PI/I80 ! converts degrees to radians 
ENDIF 
RESTORE . 

DATA 0,0,30,0,30,10,10,10,10,30,0,30,0,0,-9933,0 
COLOR 2 

read_and_plot(choice*) 

FOR i=21 _ T0 24 
PRINT AT(i,i);SPACE*(78) 

NEXT i 
LOOP 

PROCEDURE aapping_function 

• window is t-2*Yc,2*Yc] x [-Yc,Yd a square would look approx 
square 

LET x=(x+2*yc)/(4*yc 
LET y=(y+yc)/(2*yc)*(-yl)+yl 
RETURN 

PROCEDURE read_and_plot(a*) 

’ READ AND PLOT THE DATA--Adapt this to the figure you want to 
graph 
READ x,y 

translation_foraulas(choice*) 

aapping_function 

DO 

PLOT x,y 

READ x,y 

EXIT IF x=-9999 

translation_foraulas(a*) 

aapping_function 
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DRAW TO x,y 
LOOP 
RETURN 

PROCEDURE translation_forsmlas(a$) 

IF a$="M” THEN 
LET x=x+h 
LET y=y+k 
ENDIF 

IF al= B S" THEN 
LET x=x*sx 
LET y=y*sy 
ENDIF 

IF a$="R" THEN 
LET st=SIN(theta) 

LET ct=CQS(theta) 

LET xx=x 
LET x=x*ct-y*st 
LET y=xx*st+y*ct 
ENDIF 
RETURN 

tmtt*mm}m**if****mitt*»**it*i*itm 

' Another deso of using transforsations 

’ Using Translation, Scaling, Rotation subroutines 

' Using clipping routine 

' The figure progressively changes with each succeeding transforaation 

f ‘ 

view_window 

draw_axis 

i 

» LOAD UP ARRAY WITH THE PICTURE 
DIN xC100),y(100) ! Bigger than necessary 
LET i=l 

DATA 6,0,30,0,30, ljO, 10,10,10,30,0,20,0,0,-1111,0,30,20,60,30,45,50,30,30, 
DO 

READ x(i),y(i) ' " 

EXIT IF x(i)=-9999 
INC i 
LOOP 

t 

' GRAPH OUT ORIGINAL FIGURE 
LET h=o 
LET k=0 

plot figure("N“) 

COLOR 2 
DO 

• TRANSLATION MENU 

dear_textwindow 

REPEAT 

PRINT AT(i,21)|"Choose Hove Scale Rotate Suit (M,S,R 
or Q)»"j 

LET choice$=UPPER$(INPUT$(1)) 

UNTIL INSTR("NSRQ B ,choice$)>0 
PRINT choice! 

EXIT IF choice!="Q" 

IF choice$=°N” THEN 

INPUT "Enter your translation-- H,K (return) »",h,k 
ENDIF 


-3999,0 
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IF choice$=°S" THEN 

INPUT "enter the scale factors-- Sx,Sy (return) »",sx,sy 
ENDIF 

IF choice$=°R" THEN 

INPUT "Enter your rotation (degrees)—Theta (return) »",theta 
LET theta=theta*PI/l80 
ENDIF 

plot_figure(choice$) — 

LOOP 

PROCEDURE clear_textvindov 
FOR i=21 TO 24 
PRINT AT(1,i)|SPACES(70) 

NEXT i 

PRINT AT(1,21); 

RETURN 

PROCEDURE view_window 
clearjextwindow 

INPUT""Enter viewport Xs,Xl,Ys,Yl » ",xs,xl,ys,yl 

INPUT "Enter window A,B,C,D (Try -100,100,-100,100)>> B ,a,b,c f d 

COLOR 1 

BOX xs,ys,xl,yl 
RETURN 

PROCEDURE draw.axis 
IF a<0 AND b>0 THEN 
screen_coordinates(0,c) 

PLOT xv,yv 

screen_coordinates(0,d) 

BRAN TO xv,yv 
ENDIF 

IF c<0 AND d>0 THEN 
screen_coordinates(a,0) 

PLOT xv,yv 

screen_coordinates(b,0) 

DRAM TO xv,yv 
ENDIF 
RETURN 

PROCEDURE screen_coordinates(x,y) 

LET xv=(x-a)/(b-a)*(xl-xs)+xs 
LET yv=(y-c)/(d-c)*(ys-yl)+yl 
RETURN 

PROCEDURE plot_figure(a$) 

LET i=l 
DO 

EXIT IF x(i)=-9999 
IF x(i)=-llll THEN 
INC i 
ENDIF 

transfora_point(x(i),y(i),a$) 

LET x(i)=xt 
LET y(i)=yt 

screen_coordinates(x(i),y(i)) 

LET xl=xv 
LET yl=yv 
INC i 
DO 



EXIT IF x(i)=-lill OR x(i)=-9999 
transfors_point(x(i),y(i),a$) 

LET x(i)=xt 
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LET y(i)=yt 

screen_coordinates(x(i), y(i)) 

LET x2=xv 
LET y2=yv 
clip_line 
LET xl=x2 
LET yl=y2 
INC i 
LOOP 
LOOP 
RETURN 

PROCEDURE transfora_point(x,y,a$) 

IF a$= s r THEN 
xt=x+h 
yt=y+k 
ENDIF 

IF a$="S“ THEN 
LET xt=x*sx 
LET yt=y*sy 
ENDIF 

IF a$=“R" THEN 
LET st=SIN(theta) 

LET ct=COS(theta) 

LET xt=x*ct-y*st 
LET yt=x*st+y*ct 
ENDIF 
RETURN 

' serge in the clipping subroutine here 
’ PROCEDURE dip.line 

mj***j»**«i*i«**m«m*»»«***i*i**i* 

’ Snowdeao 

1 Another deao of using transforaations 
' Using Translation, Scaling, Rotation subroutines . 

’ Using clipping routine 

' One leg of a snowflake is drawn, rest is drawn using 
’ transforaations. 

' Hill use the saae prograa outline as in previous prograa. 
DIM x(13),y(13) 

DO 

' the data will be created “randomly" 

t 

view_window 

» 

’ LOAD UP ARRAY HITH THE PICTURE 
FOR i=l TO 6 
LET x(i)=i-l 
LET y(i)=2*RND-l 
NEXT i 
LET x(7)=6 
LET y(7)=0 

’ reflect for rest of snowflake arm 
FOR i=6 DOWNTO 1 
LET sx=l 
LET sy=-l 

transfora_point(x(i),y(i),"S") 

LET x(14-i)=xt 
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LET y(14-i)=yt 
NEXT i 
» 

COLOR 2 

FOR angle=0 TO 300 STEP 60 
LET theta=angle*PI/180 
transforB_point(x(l),x(2),"R") 
screen_coordinates(xt,yt) ” 

LET xl=xv 
LET yl=yv 
FOR i=2 TO 13 

transfors_point(x(i),y(i),°R B ) 
screen coordinates(xt,yt) 

LET x2=xv 
LET y2=yv 
dip_line 
LET xl=x2 
LET yt=y2 
NEXT i 
NEXT angle 

LOOP 

PROCEDURE dear_textwindow 
FOR i=2i TO 24 
PRINT AT(l,i)jSPACE$(70) 

NEXT i - . - 

PRINT AT(1,21); ' .. 

RETURN . : 

PROCEDURE view_window 
clear_texlwindow 

INPUT "Enter viewport Xs,Xl,Ys,Yl » a ,xS|Xl,ys,yl 
COLOR 1 

BOX xs,y5,xl,yl 

INPUT "Enter window - A,B,C,D(Try -10,10,-10,10) » 

RETURN . 

PROCEDURE screen_roordinatesTx,yT- 
LET xv=(x-a)/(b-a)*(xl-xs)+xs 
LET yv=(y-c)/(d-c)*(ys-yl)+y1 

RETURN 

PROCEDURE transfora_point(x,y,a$) 

IF a$="H" THEN 
xt=x+h 
yt=y+k 
ENDIF 

IF a$="S" THEN 
LET xt=x*sx 
LET yt=y*sy 
ENDIF 

IF a$="R" THEN 
LET st=SIN(theta) 

LET ct=COS(theta) 

LET xt=x*ct-y*st 
LET yt=x*st+y*ct 
ENDIF 

RETURN 

PROCEDURE dipjine 

» MERGE IN YOUR LINE CLIPPING SUBROUTINE 
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Program Seahavk 
' Seahawk 
' by Dick Plagge 
' Highline CC 

' A PR08RAH THAT USES THE TRANSFORMATION MATRIX 

f 

DATA 0,0,15.3,0,15.2,1,15.2,2,15.5,3,15.8,4,16,4.4,18.1,4.7 
DATA 14,8,0,8,0,0,-1111,0,0,10,14,10715,10,15.5,11,16,11.3 
DATA 17,12.1,18,12.7,13,13.2,20,13.6,21,14,22,14.1,23,14.2 
DATA 23.5,14.22,24,14.2,25,14,26,13.7,27,13.3,28,12.8,29,12.3 
DATA 29.5,12.25,30,12.3,30.5,12.6,31,13,31.6,13.8,31.8,13.5 
DATA 31.5,13,31,12.4,30,12,29,11.85,28,12,27,12.4,26,12.8,25 
DATA 13.1,24,13.25,23,13.3,22.5,13.33,22,13.3,21,13.2,20,13,19 
DATA 12.6,18.5,12.2,18,11.8,17.5,11.4,17,10.8,16.5,5.5,17,4.8 
DATA 16.5,3.9,16.2,3,16,2,16,1,16.1,0,39,0,40,-.1,41,-.3,42.2 
DATA -.6,42.7,.4,42,.7,41,1,39,1.6,38,1.85,37,2.2,36,2.3. 

DATA 35,2.5,34,2.6,32,2.8,28,3.1,25,3.1,24,3,23,2.8,22,2.5 
DATA 21,2,20.5,1.4,20.2,1,19.7,1,19.8,1.5,20,2,20.5,2.6 
DATA 21,3.2,21.5,3.7,22,4,23,4.2,24,4.3,25,4.4,27,4.5,29 
DATA 4.3,32,4,35,3.5,37,3,40,2,42.8,.9,43,1.25,43.2,2,43.6 
DATA 4,43.8,5,43.9,5,43.7,8,43.5,9,42.9,10,43,11,41,11.5 
DATA 40,11.8,33,12,35,12,34,10,33.2,8,33,7.7 
DATA 32.6,7,32,6.3,31.5,5.85,31,5.7,30,5.5,28,5.6,26,5.8 
DATA 25,6,24,6,23,5.8,22,5.55,21,5.2,20,4.5,13,3.5,18.5 
DATA 2.5,18.2,2,18,1,19,1,19.4,2,20,3,20.5,3.6,21,4,22,4.5 
DATA 23,4.6,26,5,28,4.3,30,4.7,33,4.2,35,3.9,38,3,41,2,43,1.25. 
DATA 42.8,.9,42.7,.4,42.2,-.6,41,-3.7,42,-3,43,-2.2,44,-1,45 
DATA .5,46,3,46.5,4.5,46.8,6,46.9,8,46.8,9,46.6,10,46.4,11 
DATA 46,12,45,13,44,13.6,43,14,24,14.2,35,14.2,33,16,32,16.6 
DATA 30,17.6,29,17.8,28,18,0,18,0,10,-1111,0 
DATA 30,10,28,11.3,26,12,24,12.4,23,12.6,22,12.6,21,12.5,20 
DATA 12.2,19,11.7,18,11,17.3,10,17,9.3,16.3,9,16.8,8,16.9 
DATA 6.5,17,6,17.5,5.4,18,5,18.5,6.3,19,6.5,20,5.8,21.5,7,22.5 
DATA 6.9,23,6.8,24,5.65,25,6.6,26,6.65,27,6.8,27.5,7,28,7;3 : ’ 

DATA 28.7,8,29.4,9,30,10,-1111,0 

DATA 28,10,27,10,26,10.4,25,10.9,24,11.3,23,11.4,22.5,11.3,22,11.2 
DATA 21.5,11,21,10.7,20,10,18.4,10,18.5,9.5,20,9.5,20.5,9.3 
DATA 21,10.3,22,10.8,23,11,24,10.9,25,10.3,26,10,27,9.6,28,9.6 
DATA 28,10,-1111,0,28,9,28,9.4,27,9.3,26,8.8,25,8.2,24,7.8 
DATA 23.5,7.7,23,7.8,22,8,21,8.5,20.5,8.8,20,9,28.5,9,18.5,8.6 
DATA 20,8.6,20.5,8.5,21,8.3,21.5,8,22,7.7,23,7.5,23.5,7.4,24,7.5 
DATA 24.5,7.6,25,7.8,25.5,8.2,26,8.5,26.5,8.7,27,8.8,27.5,8.9,28,9 
DATA -1111,0,21,9,22,9,22,9.3,21,9.3,21,9,-1111,0 
DATA25,9.5,24.95,9.9,24.8,10.3,24.6,10.6,24.3,10.8,23.9,10.95,23.5 
DATA 11,23.1,10.95,22.8,10.8,22.4,10.6,22.2,10.3,22.05,9.9,22,9.5 
DATA 22.05,9.1,22.2,8.7,22.4,8.4,22.8,8.2,23.1,8.05,23.5,8 
DATA 23.9,8.05,24.3,8.2,24.6,8.4,24.8,8.7,24.95,9.1,24,9.5 
DATA -9999,0 
REPEAT 

FOR i=21 TO 24 
PRINT AT(l,i);SPACE$(70); 

NEXT i 

PRINT AT(1,21); 

PRINT “Enter entries of satrix I A C 
PRINT “ I B D 


INPUT “Enter A,8,C,D- 
COLOR 1 
DEFFILL 0,1 


try 5,0,0,2 as a start ",a,b,c,d 
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PBQX 0,0,639,159 
COLOR 2 
RESTORE 
DO 

READ x,y 
natrixjrans 
PLOT xjy 
DO 

READ x,y 

EXIT IF x=-llll OR x=-9999 
GOSUB aatrix_trans 
DRAM TO x,y" 

LOOP 

EXIT IF x=-9999 
LOOP 

UNTIL a=0 AND b=0 AND c=0 AND d=0 

i 

PROCEDURE aatrix trans 
LET xl=x 
LET yl=y 
LET x=a*xi+b*yl 
LET y=c*xl+d*yl 
LET x=x+200 
LET y=-l*y+100 

RETURN 

************************************************* 


sl programing a s sj 1 gnme nt 

Modify the Seahawk program so you enter in the matrix that 

translation^ie. Homogeneous Coordinates ). 


allows 
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chapter ten 

SURFACES OF REVOLUTION 

This chapter will ease us into the world of 3-dimension. A very 
striking figure is_ a three dimensional graph that is a symmetric 
about the origin. The figure below is a sample of such a graph. 

I 



The figure is formed by considering the graph of a simple 
function in the x-y plane and then rotating that figure around the 
y-axis.(See the sketch below) This is called a surface of 
revolution. 



We will imagine a coordinate system in which the z-axis comes 
out of the screen toward the viewer. (We will discuss coordinate 
systems in much more detail in the next chapter.) I will draw the 
surface in cross-sections parallel to the x-y plane(screen) from 
back to front. Imagine we're looking down on the surface from the 
top.(See the two figures on the next page.) If R is the radius of 
revolution, then the z-values will go from -R to R. 
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Program to graph surface of revolution: 

1 ThreeD 

' A surface of revolution prograe 
i I 

LET x1=639 
LET yl=199 
LET xc=INT((xl+l)/2) 

LET yc=INT((yl+i)/2) 

LET r=xc-5 ! sake radius of rotation slightly saaller than 1/2 screen 
FOR z=-r TO r STEP 10 ! picks cross-sections froa back to front 
LET xffl=S8R(r*r-z*z) 

FOR x=-xa TO x* i picks x values for each cross-section 
LET t=SSR(x*x+z*z) ! value to put into function 
LET y=10*C0S(0.05*t)+300/(t*t+3) ! THE FUNCTION HE'RE PLOTTING 
* HERE NOT USING A FANCY HINDOHING FUNCTION, THE VALUES OF 
1 X AND Y ARE BASICALLY SCREEN COORDINATES 
LET y=y-0.09*z ! tips back edge up—corresponds to about 5 degrees 
LET xl=x+xc 
LET y=-i*y+yc 
IF y<yl AND y>0 THEN 
PLOT xl,y 
COLOR 0 

DRAH xl,y+l TO xl,yl 
. COLOR 1 
ENDIF 
NEXT x 
NEXT z 



The value of xm is the maximum x 
for that cross-section. The first 
top view shows the relation between 
z, r, xm. By the Pythagorean 
theorem, xm=SQR(r a -z a ). 

The second top view shows an x 
picked on the cross-section. The 
distance t to that point (from the 
origin) is related to z and x 
Again the Pythagorean theorem comes 
into action. This time as 

t = SQR(x z +Z 52 ) . 

Any point a distance of t from the center of rotation will have 
the same y-value, so consider t the value to be placed into the 
function, and compute f(t). 
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The sneaky part of the- program is the way we do not show any 
hidden lines. Since we're drawing the surface from back to front, 
any cross-section in front will automatically block out "stuff" 
behind it. So the program plots the point, changes the color to 
the background color, then draws a line from below that point to 
the bottom of the screen. This effectively erases the previously 
drawn "stuff" that was lower than the plotted point. The color is 
then changed back to the plot color. 

We now tip the back of the figure up and the front down for 
better viewing. If z is negative, then y-0.09x is greater than 
the original y since -.09*z would be positive. Similarly if z is 
positive, y-0.09z makes y smaller. The 0.09 gives about a 5 
degree tip. 

ANOTHER WAY . .. . 

A second approach to graphing the same surface is to graph it 
from front to back. At each x-postion keep track of the larges t 
y-value plotted and the smallest v-value plotted. Use two arrays 
(ub and lb—DIM. ub(xl ), lb(xl) ) to keep track of this information. 
If a point is above or below a previously drawn point then plot 
it. See the next program for the implementation of this idea. 
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’ ThreeDee 

’ Another surface of revolution prograa 


LET xl=639 
LET yl=199 
DID ub(xl),lb(xl) 

ARRAYFILL ub(),-9999 
ARRAYFILL lbO,9999 
LET xc=INT((xl+l)/2) 

LET yc=INT((yl+l)/2) 

LET r=xc-5 ! aake radius of rotation slightly saaller than 1/2 screen 
FOR z=r TO -r STEP -10 ! picks cross-sections froa front to back 
LET xa=INT(SQR(r*r-z*z)) 

FOR x=-xa TO xa ! picks x values for each cross-section 
■ LET t=SQR(x*x+z*z) ! value to put into function 
LET y=lO*CQS(0.05*t)+300/(t *t+3) ! THE FUNCTION HE'RE PLOTTING 
» HERE NOT USING A FANCY HINDOHING FUNCTION, THE VALUES OF 
' X AND Y ARE BASICALLY SCREEN COORDINATES 
LET y=y-0.Q3*z ! tips back edge up—corresponds to about 5 degrees 
LET xl=x+xc 


LET y=-l*y+yc 
LET upper=ub(r+x) 

LET lower=lb(r+x) 

IF y<=upper AND y>=lower THEN 
! do nothing 
ELSE 

ub(r+x)=MAX(upper,y) 
lb(r+x)=HIN(lower,y) 

IF y<yl AND y>0 THEN 
PLOT xl,y 
ENDIF 
ENDIF 
NEXT x 
NEXT z . 

"SGET picture* isaves screen in the variable picture! 

PRINT AT(1,22); 

INPUT "Do you wish to print this picture? “,answer$ 

IF UPPER*(answer*)=“Y“ THEN 

PRINT AT(1,24)j“Be sure printer is on. Press ALTERNATE HELP- 
PRINT “ when this aessage disappears. Hhen done, press RETURN 0 
PAUSE 250 


SPUT picture* !puts picture back on screen 
INFUT “”,a* ! keeps prograa froa ending 
ENDIF 

INPUT “Do you wish to save this picture? ",answer* 

IF UPPER*(answer*)=“Y“ THEN 
INPUT “Enter nase of file to save this under. “,naae* 

SPUT picture* ! put picture on screenfno writing!) 

BSAVE naae$,XBI0S(2),32000 ! Xbios(2) gives address of screen 
» It takes 32000 bytes for a picture 


ENDIF 


' The following shows how to put picture on screen froa the file 
’ BL0AD naae*,XBI0S(2) ! you eight want an input to get the 
' naae of the file 
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The second method seems to be faster although both are very 
time-consuming. The first method uses a lot less memory since it 
is "expensive" byte-wise to store numbers in arrays. But what the 
heck—we got memory to burn these days! 

Try other functions in the program!! You might have to 
experiment a bit to get the constants adjusted so you get a figure 
you like. You might consider increasing the step size on the 
outside loop as your experimenting, so you can get a rough idea 
quickly as to how your creation looks. 

By the way, the variables "upper" and "lower" are upper and 
lower bounds of the values of y at each x-position. Since the 
screen coordinates are somewhat upside down, these are the reverse 
of their physical location on the screen. 

Later, after we've discussed the process of graphing in 3-D more 
thorougly. I'll try to present a fairly simple program to graph 
general 3-D surfaces. 
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chapter eleven 

"THE REAL" REAL WORLD 
OF 3 —D 

The job at hand(and for the rest of this book) is to look at the 
three dimensional world. We need to figure out how to represent 
3-D objects on the little two dimensional screen of the computer. 
Making the 3-D object look realistic can involve such diverse 
topics such as perspective, shading, color, hidden lines and 
surfaces, texture, and shadows. At this point, only the most 
sophisticated computers have the capability to represent all of 
these characteristics realistically. This book will stick to the 
most basic of 3-D graphics. Of course, our "old" topics of 
windows, viewports, and clipping are still relevant. 

The mathematics needed to discuss all of these topics is much 
beyond a short book in graphics. Also, the equipment needed to 
demonstrate all of these topics is beyond the resources of the 
average person with a microcomputer.(However, the advent of the 
new cheap 16-bit computers with "super" high resolution graphics 
is putting this into reach.) 

To start off, we need to be aware of how points in 3-D can be 
represented. Some of the standard coordinate systems will now be 
looked at. 
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As you can see, our 3-D rectangular coordinate system is a 
natural extension of the 2-D cartesian coordinate system. Many 
formulas are just extensions of 2-D formulas. For example, the 
distance between two points (Xi,yi) and (x=:,ya:) is given by 
/(x 2 -xi) 2 +(y 2 -yi) 2 . In three dimensions, the distance 

between (xi ,yi ,Zi) and (x 2 ,y 2 ,z = ) is given .... by 

/(X 2 -X 1 ) 2 +(y^-yi) 2 +(Z 2 -Z 1) 2 . The orientation of the xyz-axis 

can vary from book to book, however the one shown above is the 




most common in mathematics books. 

However the one to the right is a common 
orientation for computer graphics. Note 
that the x-y axis are in normal position 
as in the two dimension case, but the z 
axis comes out toward us.(Out from the 
screen, if the screen represents the x-y 
plane.) ^2 

If you can imagine a stick model of the axis, other 
representations can be imagined by simply rotating it into 
different orientations. All are what called "right-handed" 
systems. The name stems from using the fingers of the right hand 
to represent the x-y axis and the thumb to represent the z-axis. 
There are a couple of ways to do this. Both, however, strain my 
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artistic abilities to draw pictures. 

or 

Imagine pointing all your fingers 
toward the x-axis, then curl your fingers 
into a fist toward the y-axis. In doing 
this, your thumb will point in the 
direction of the z-axis. 

In a left-handed system the z-axis would point in the opposite 
direction. It is not entirely uncommon to see the z-axis point 
into the screen in some graphics books. This is often the case 
when you wish z to represent the concept of depth. 

Also, a convention has to be established about what is a 
positive angle when considering rotation. 



The convention we shall use is to look 
the axis you're rotating about, toward 
origin. Rotation counter-clockwise 

considered positive. 



-> > 

* ^ 


CYLINDRICAL COORDINATE SYSTEM 


a £ 




The cylindrical system is basically 
polar coordinates with the z component 
added for the 3rd dimension of 
height.(The point P' is a point in the 
x-y plane.) The name "cylindrical 
coordinate system" comes from the fact 
it is very easy to express the 
equation of a cylinder in this system. 


Consider the set { (r,0 ,z): r=2 } 

Imagine a point P in that set, where 0 can be any angle, z can 
be any height, but the point always has to be 2 units away from 
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the z-axis. If you visualize an "infinitely" long cylindrical 
tube centered around the z-axis, you've got it!! This undoubtedly 
is the basis of the name of this system. 

SPHERICAL COORDINATE SYSTEM 


In actuality, the spherical 
coordinate system is more (in 
spirit) like the 2-D polar 
coordinate system. As you can see 
from the diagram of the system, we 
measure the distance to the point 
directly from the origin. 

Angles are then used to describe the orientation of that line 
segment. Since we're in three dimensions we have an extra degree 
of freedom, so it takes two angles to fix its location. By 
convention, $ is an angle such that 0<C f> tc and p >0. 

The name "spherical" comes from the fact that a sphere can 
easily be described by { (p,9 ,f> ): p =10 } for example. 

.Conversion from cylindrical or spherical coordinates to 
rectangular coordinates is often a necessity. (It may be 
theoretically more elegant to use cylindrical or spherical 
coordinates to describe a situation, but if it has to be adapted 
to a computer system that "knows" only rectangular coordinates, 
then conversion is necessary.) Only elementary trig relationships 
are necessary to make the conversions. The following list come 
easily by considering the different right triangles involved. 


*2 



P(-?>,<f> 


I 

I 




Cylindrical 


Spherical 


(r,0,z)—>(x,y,z) 


(p,9,$)—>(x,y,z) 


x = rcosG 
y = rsin9 
z = z 


x = psin$cos0 
y = psinf>sin0 
z = pcosf> 
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pr oblems 

1 . 

(a) Establish the formula for calculating distance of a point 
from the origin in rectangular coordinates. 

(b) Find the distance from the origin of P:(2,4,10) 

(c) What is the distance between P:(2,4,10) and Q:(-2,7,9)? 

•2. Draw a couple of other orientations of right handed coordinate 
systems other than what is shown in the text. 

3. Find the rectangular coordinates of the following points: 

(a) (5,it/6,10) in cylindrical coordinates. 

(b) (7,7 tc/ 4,-2) in cylindrical coordinates. 

(c) (7 ,tc/6,tc/ 4) in spherical coordinates. 

(d) (10,5 tc/6,tc/ 3) in spherical coordinates. 

4. Describe -in words -the figures these equations describe. Try 
to sketch a picture of each. 

(a) { (r,9,z) :z=5 } (b) { tp,9,4>): 4> = it/4 } 

(c) { (-f,9,4> ) :^cos$ =4 } (d) { (r,9,z) : 9 = irr/4 } 
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chapter ■twelve 

DISPLAYING 3—D 
ON A 2—D SCREEN 



I'll gradually work into this topic. I'll make some simplifing 
assumptions first, look at a "easy” way to represent 3-D objects 
then work up to a more realistic representation. 

Now, normally, when coordinates (x,y,z) are given(especially in a 
math book) it is with respect to the following d 

However, our screen is 
usually thought of as: 

Whereas the z-coordinate is often used to represent vertical 
height of a point above the xy plane, we would normally use y to 
represent height in a two dimensional situation. So, to simplfy 
things in my programs I will think of 3-D as in the diagram below. 

I will think of the TV -, screen - 
being positioned on the positive z- 
axis and the viewer being further 
out on the z-axis. The TV screen 
will be perpendicular to the z-axis 
and parallel to the xy plane. 

To further simplify things it will be assumed the z-axis will 
penetrate the sreen in the center. 

The viewport will be the entire screen and the window is 
basically the same as the viewport, but with the origin translated 
to the center (ie. window = C-xl/2,xl/2]x[-yl/2,yl/2] viewport = 
t 0,xl ]x[ 0,yl ] ). The coordinates of the 3-D objects will be 
chosen so that clipping will not normally be a problem. 
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Let's write programs to display a 



PARALLEL ORTHOGRAPHIC PROJECTION 


"solid" stairstep object: 

A typical vertex of the 
object will have coordinates 
(x,y,z), the "depth being given 
last. 


The simplest type of representation of a 3-D 
object is the parallel orthographic 
projection. (Say that 5 times fast!!!) This 

type of projection simply ignores the depth . I-, 

information. (ie. (20,0,10) becomes (20,0); 

(20,0,0) becomes (20,0) also; etc.) ‘ The 
stairstep would look like the diagram to the 
right (looking down the z-axis ). 

If we put our'TV screen on the x-axis( -and. 
made the x coordinate represent depth, we 
would get the following view: (the x 
coordinate would be ignored in graphing the 

point, and only the y and z information would 1 - 1 

be relevant) 

These are views you would see on a blueprint in certain drafting 
applications. Given several "side" views one can learn to 
interpret what the 3-D object must look like. 

See the program at the end of this chapter demonstrating 
parallel orthographic projections. It will give you a "front" 
(looking down the z-axis) view, a "side" (looking down the x-axis) 
view, and a "top" (looking down the y-axis) view of the object. 
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You could add more objects at the end of the program. Many 
variations to this program could be added for your own purposes 
and pleasure. 

The subroutines used are fairly obvious. The subroutine called 
Position_object places the object on the screen. The H and K 
allow you to position various views at different places on the 
screen. It is a combination windowing function and translation 
routine. The subroutine called change_views "renames" the 
coordinate variables for the various views.(Remember, everything 
is "x and y" to the computer screen.) 

Can you determine what objectl is before you run the second 
program? You may find it hard to determine an objects shape from 
these three parallel projections, but I'm told that with practice 
it's possible! 

PERSPECTIVE PROJECTION 

It's pretty obvious that parallel orthographic projections are 
not what most people consider as a 3-D representation of an 
object, however. Generally, one usually wants to include 
perspective in representing an object to give an illusion of 
realistic depth. The idea of perspective is that as one moves an 
object further away from the viewer, it will appear to be smaller. 
Thus the back face of a cube will appear smaller than the front 
face, railroad tracks will appear to be closer together as one 
views them receeding, etc. There are different types of 
perspective schemes, but I won't bother to identify all the 
different types in this book. 
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To see how a point projects on the screen(under our simplifing 
assumptions made earlier) consider the following diagrams: 
( diagrams made by "god” observing the object, the axis, and the 
observer) 




Considering similar triangles in the above pictures: 

x' x x(de-ds) y' y y(de-ds) 


de-ds de-z de—z de-ds de—z de—z 

As you can see from the diagram, the distance to the eye(De) and 
the distance to the screen(Ds) are measured from the origin of our 
three dimensional space( from which-our object is referenced ). 
This is not always the most natural from the viewer's point of 
view, but in the interest of the famous kiss principle I'll not 
try to interact a viewer's coordinate system with the object's 
coordinate system. 

Now examine the second program at the end of this chapter. I 
tried to keep it almost the same as the previous program. The 
subroutine perspective_window contains the perspective formulas 
developed previously and incorporates the windowing function. The 
subroutine scale_translate_object scales then translates the 
object before we view it. 

You have to play with the values for de and ds to get the 
perspective you want. Also be careful of the scale factor or 
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you'll find yourself inside the object. Also you might find some 
points being off screen. This won't really be a problem for this 
computer, as GEM usually will "plot points off screen" with . no 
error messages. This is not true universally!! If you modified 
this program to work within a smaller viewport other than the 
whole screen, then you would want to think about incorporating the 
clipping subroutine. 

Since both of my objects built into the program are "tied" to 
the origin, I've applied the scaling routine first, then 
translation. I think that will seem to give the desired results 
most naturally. 

There are no hidden lines, so sometimes the lines in "back" 
confuse the eye. Sometimes you really have to blink the eyes to 
make it look the way you know it.has to look! 

Feel free to modify the programs, put in your own objects. 
Other books may help you to jazz up the programs so that hidden 
lines and faces can be determined and not drawn. The hidden line- 
face problem is not a simple.one, . and I've decided not to .tackle 
it in this introductory "book. 
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' A siaple deao of three-diaensional objects 
’ Parallel Orthographic Projection 

objectO: - - . 

DATA 0,0,10,20,0,10,20,20,10,10,20,10,10,30,10,0,30,10,0,0,10,-1111,0,0 
DATA 0,0,0,20,0,0,20,20,0,10,20,0,10,30,0,0,30,0,0,0,0,-1111,0,0 
DATA 20,0,10,20,0,0,-1111,0,0 

DATA 10,20,10,10,20,0,-1111,0,0,10,30,10,10,30,0,-1111,0,0 
DATA 0,30,10,0,30,0,-1111,0,0,0,0)10,0,0,0,-9999,0,0 
objectl: 

DATA 0,0,0,0,0,50,50,0,50,50,0,0,0,0,0,25,25,25,0,0,50,-1111,0,0 
DATA 50,0,0,25,25,25,50,0,50 
DATA -9999,0,0 

r 

LET xc=320 

LET yc=80 ! Viewport [0,2*xclx[0,2*yc] 

LET b=200 

LET d=100 ! window [-b,blx[-d,d] 

f 

DO 

CLS 

draw_axis 

set_up_textwindow 

INPUT "enter object # (0 or 1) to graph(-99 to end) »",object 
EXIT IF object=-99 

- i 

set_up_te'xtwindow 

PRINT "The front view, no translation, no scaling" 

COLOR 1 

plot_object(l,l,l,0,0,"F") 
firsttiae!=TRUE 
PAUSE 200 
DO 

set_up_textwindow 

‘INPUT "EnterScale factor—Sx f Sy,Sz »",sx,sy,sz 
INPUT "Enter-Translation-horizontal,vertical »",h,k 
INPUT "Enter view—Front,Side,or Top(F,S,or T) »°,view$ 

IF firsttiae! THEN 
COLOR 0 ! Erase original figure 
plot_object(l,l,l,0,0,"F") 
draw axis 
firsttiae!=FALSE 
ENDIF 
COLOR 1 

plot_object(sx,sy,sz,h,k,UPPER$(view$)) 
set up textwindow 

PRINT "CHANGE OBJECTS OR CLEAR SCREEN?(y,or n) "j 
answer$=INPUT$(l) 

EXIT IF UPPER$(answer$)="Y" 

LOOP 
LOOP 

PROCEDURE draw.axis 
COLOR 1 

DRAW 0,yc TO 2*xc,yc i x axis 
DRAW xc,0 TO xc,2*yc ! y axis 
PRINT AT(77,ll);b 
PRINT AT(40,l)jd 
RETURN 
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PROCEDURE restore.data 
IF objects THEN 
RESTORE objectO 
ELSE 

RESTORE object 1 
ENDIF 
RETURN 

PROCEDURE 5 et_up_textwindow 
FOR i=21 T0’25" 

PRINT AT(l,i)jSPACE$(79)j 
NEXT i 

PRINT ATCl,21)j 
RETURN 

PROCEDURE plot_object(x5cale,yscale,zscale,horiz,vert,view$) 
restore_data 
DO 

READ x,y,z 

scale_object(x5cale,yscale,zscale) 

change_views(view$) 

position object(horiz,vert) 

PLOT x,y” 

DO 

READ x,y,z 

EXIT IF x=-9999 OR x=-llil 
scale_object(xscale,yscaIe,zscale) 
change_views(view$) 
position_object(horiz,vert) 

DRAW TO’x.y 
LOOP 

EXIT IF x=-9999 
LOOP 
RETURN 

PROCEDURE scale_object(x5cale,yscale,z5cale) 

LET x=x*xscale . ' 

LET y=y*yscale 
LET z=z*zscale 
RETURN 

PROCEDURE change_views(vie«$) 

IF vieu$=°S" THEN 
LET x=z 
ENDIF 

IF viewi^T" THEN 
LET y=z 
ENDIF 
RETURN 

PROCEDURE position_object(horiz,vert) 

LET x=x+horiz ! Apply the translation 
LET y=y+vert 

LET x=(x+b)/b*xc ! Hindowing function 
LET y=(d-y)/d*yc 
RETURN 
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' A siaple demo of three-diaensional objects 

' With Perspective 

objectO: 

DATA 0,0,10,20,0,10,20,20,10,10,20,10,10,30,10,0,30,10,0,0,10,-1111,0,0 
DATA 0,0,0,20,0,0,20,20,0,10,20,0,10,30,0,0,30,0,0,0,0,-1111,0,0 
DATA 20,0,10,20,0,0,-1111,0,0 

DATA 10,20,10,10,20,0,-1111,0,0,10,30,10,10,30,0,-1111,0,0 
DATA 0,30,10,0,30,0,-1111,0,0,0,0710,0,0,0,-9939,0,0 
objectl: 

DATA 0,0,0,0,0,50,50,0,50,50,0,0,0,0,0,25,25,25,0,0,50,-1111,0,0 
DATA 50,0,0,25,25,25,50,0,50 

DATA -9339,0,0 
» 

LET xc=320 

LET yc=80 ! Viewport [0,2*Xc]x[0,2*YcI 
LET b=200 

LET d=100 ! window [-B.BlxE-D,D3 

i 

DO 

CLS 

draw_axis 

set_up_textwindow 

INPUT “enter object I (0 or 1) to graph(-99 to end) »“,object 
EXIT IF object=-93 

i 

set_up_textwindow 

LET~de=300 . - 

LET ds=100 

PRINT “Scale of 1, no translation, distance to eye is °;de 
PRINT “Distance to screen is “jds 
COLOR 1 

plot_object(l,0,0,0) 
firsttime!=TRUE 

PAUSE 200 - - 

DO 

set_up_textwindow 

INPUT “Enter overall Scale factor--Sf »“,sf 
INPUT “Enter Translation in X and Y and Z direction »“,h,k,l 
PRINT “Present distance to eye, distance to screen is °;de,ds 
INPUT "Enter distance to eye,distance to screen » "jde,ds 
IF firsttime! THEN 
COLOR 0 ! Erase original figure 
plot_object(l,0,0,0) 
draw axis 
firsttime!=FALSE 
ENDIF 
COLOR 1 

plot_object(sf,h,k,l) 

set_up_textwindow 

PRINT “CHANGE OBJECTS OR CLEAR SCREEN?(y,or n) "} 
answer$=INPUT*Cl) 

EXIT IF UPPER$(answer$)="Y“ 

LOOP 

LOOP 

PROCEDURE draw.axis 
COLOR 1 

DRAW 0,yc TO 2*xc,yc ! x axis 
DRAW xc,0 TO xc,2*yc ! y axis 
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PRINT AT(77,lt)jb 
PRINT AT(40jl)jd 
RETURN 

PROCEDURE restore.data 
IF objects THEN 
RESTORE objectO 
ELSE 

RESTORE objectl ~ 

ENDIF 

RETURN 

PROCEDURE set up_textuindow 
FOR i=21 T0’25* 

PRINT AT(l,i);SPACE$(79)j 
NEXT i 

PRINT AT(1,21)j 
RETURN 

PROCEDURE pIot_object(scale,h,k,1) 
restore_data 
DO 

READ x,y,z 

scale_translate_object(scale,h,k,1) 
perspectivejjindou 
IF NOT inside figure! THEN 
PLOT x,y 
ENDIF 
DO 

READ x,y,z 

EXIT IF x=-9339 OR x=-IIli 
scale_translate_object(scale,h,k,I) 
perspective_yindow 
IF NOT inside_figure! THEN 
DRAH TO x,y 
ENDIF 
LOOP . 

EXIT IF x—-3999 
LOOP 
RETURN 

PROCEDURE scale_translate_object(scale,h,k,l) 
LET x=x*scale+h 
LET y=y*scale+k 
LET z=z*scale+l 
RETURN 

PROCEDURE parspective_uindov 
LET inside figure!=FALSE 
LET dl=de-z 
IF dl>0 THEN 
LET x=x*(de-ds)/dl 
LET y=y*(de-ds)/dl 
ELSE 

LET inside figure!=TRUE 
ENDIF 

LET x=(x+b)/b*xc ! Windowing function 
LET y=(d-y)/d*yc 
RETURN 
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EQUATION OF A LINE 

As in the two dimensional case, the parametric equations 

x = Xt+ (X 2 -Xi)t 
y = yi+ (ya-yi)t 
z = Zi + (Z 2 -Zi)t 

represent a line going through points Pi and Pa. If 0<t<l, then 
we have the line segment between Pi and Pa. 

It is not uncommon to define the concept of "adding" a point and 
a vector to get another point. If you look at the diagram above, 
you can see that starting with point Pi , adding to it the 
components of v, the coordinates of Pa emerge. 

(ie. Xi + (Xa -Xi )=Xa ,yi + (ya -yi )=ya ,Zi + (Za -Zi )=Za ) 

So, if P=(a,b,c) is a point and v is a vector with components 
<vt,va,v 3 >, then P + v = (a+vi,b+v 2 ,c+v 3 ). 

Now, the parametric equations can be written as (x,y,z) = P + tv 
where (x,y,z) is any point on the line. (By the way, this could 
have been done in the two dimensional case also. I just chose not 
to throw everything in all at one time!!) 
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problems 

1. What are the parametric equations of a line that contains the 
points (2,3,-4) and (-2,5,7)? 

2. What are the parametric equations of the line segment between 
(3,-4,6) and (4,3,-7)? 

3. A point (3,4,-2) lies on line L. A vector <-2,3,5> is 
parallel to the line. Find the parametric equations for L. 

DOT PRODUCT 

Let v = <vx,Va,Va> and u = <Ux,Ua,Ua> be two vectors. The dot 
product u»v is a number given by: 

U*V = UiVi + UstVa + UaVs. 

(The dot product could have been defined similarly for a two 
dimensional vector.) 

By the Pythagorean theorem, it is easy to see that the length of 
a vector v=<vi,v 2 ,v 3 > is given by v4i 2 + v 2 2 + v 3 s . [vj is - often 
the symbol used for the length of v.( j v[ is sometimes called the 
norm of the vector v.) 

Dot products are fun, and there are some easy properties to 
prove. Can you prove these? 

(a) v • v = | v 

(b) u*v = v»u (commutative property) 

(c) u°(v+w) = u»v + u»w (distributive property) 

(d) (tu)«v = t(u«v) = u»(tv) 

(e) |u•vj < |u|{v| 

(f) 1 v has length l.(a unit vector) 
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The dot product can be useful in expressing the angle between 

By the law of cosines, 

|u-v| 2 = ju|-+7v|» - 2[ufjv|cos0 

Now, simplifying the above equation using the rules about dot 
products, you get these equivalent equations: 

(u-v)'(u-v) = u*u + v*v -2|u||v|cos8 
u*u -2u«v + v«v = u*u + v«v -2|u||vjcos0 
-2u°v = -2ju[[v|cos8 

, 0<8<180* 

Since cos 90*= 0, we can formulate this rule: u and v are 

perpendicular whenever u*v =0. We’ll have an occasion to use 
this later on! 

problems 

1. Let v = <3,4,7> and u = <4,5,-4>. find u«v and |u| and |v|. 

2. Let v = <4,8,-2>. Find the length of v. 

3. Let u = <5,4,2>. Find u«u. Find a unit vector in the 

direction of u . 

4. What is the cosine of the angle between <2,-4,5>and <-3 2,-7>? 

5. Find the angle between <3,2,4> and <2,5,-4>. 

6. Let u = <3,5,-8>. Find a nonzero vector - that is 




perpendicular to u. 
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EQUATION OF A PLANE 

To fix a particular plane in space, all we need to know is a 
point Po= (Xo, y©,z©) on the plane and a vector n = <a,b,c> that is 
perpendicular to _the plane. (There are many planes through a 
point, but the vector n will fix it in only one direction.) 

Let P(x,y,z) be any point on the 
plane. The vector n will be 

perpendicular to any vector on the 
plane. In particular n is 

perpendicular to the vector from Po 
to Pi. 

or <a,b,c>«<x-x©,y-y©,z-z©>=0 
Finally, a( x-xo) +b(y-y©)+c( z-zo) =0 

This is considered an equation of a plane, containing the point 
Po and with a vector n perpendicular to the plane. (By the way, 
sometimes we say n is normal to the plane, rather than saying it 
is perpendicular.) 

pr oblems 

1. Find the equation of a plane containing (0,0,0) and is normal 
to <0,0,1>. 

2. Find the equation of the plane containing (2,-4,3) is is 

perpendicular to <-l,2,7>. 

3. Find the equation of the plane containing (3,5,7) and is 

parallel to the plane given by 2x-5y+3z=4. 

4. Find the equation of the plane containing (2,1,3), (-1,2,1), 

and (1,-2,2). (This is harder than the previous 3 problems!) 



So, n*<x-xo,y-yo,z-Zo> =0 
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CROSS PRODUCT OF TWO VECTORS 

Knowing how to find a vector that is perpendicular to two given 
vectors u and v is a handy thing when working in 3-D geometry. 
The vector we will_ find that is perpendicular will be called the 
cross product (denoted by u x v ) of the vectors involved. 
Let's develop a way to find a perpendicular vector. 

Let u = <ui,u 2 ,u 3 > and v = <vi,v 2/ v 3 > be two vectors. Suppose 
the vector n = <a / b / c> is to be perpendicular to u and v. For 
this to be true we would have n u=0 and n v=0. 
aui + bus + cua =0 

and 

avi + bvs 2 + cv 3 = 0 
or equivalently au 2 + buz = -cu 3 

avi + bva =.-c-Vs 

Solve these equations simultaneously for a and b: 

(mult, by Vi ) auiVi + buaVx = -cu 3 Vi 

(mult, by -Ui ) -auiVi - buiv 2 = cuiv 3 

b(U 2 Vi - UiV 2 ) = C(UiV 3 - UaVl) 

UiVa-UaVi • UaVx^UxVa 

b = - c = - —— C 

U 2 Vi-UiV 2 UiV 2 -U 2 Vi 

Similarly, solve for a. 

(mult, by v 2 ) auiv 2 + bu 2 v 2 = -cu 3 v 2 

(mult, by -u 2 ) -au 2 Vi - bu 2 v 2 = cu 2 v 3 

a(uiv2 - u 2 Vi) = c(u 2 v 3 - u 3 v 2 ) 

u 2 v 3 -u 3 v 2 
a = - c 

Uiv2-u 2 Vi 

There are infinitely many vectors perpendicular to a given pair 
of vectors(could have various lengths, be 180 degrees from each 
other). This wide variety of different vectors is expressed in 
the above equations in that there are infinitely many values for a 
and b, depending on what you choose for the value of c. It turns 
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out 

that a 

particularly good 

choice for c 

is to let c 

be the 

value 


o£ 

the 

denominators in 

the above 

fractions. 

Thus 

let 


c 

= U 1 V 3 

U2V1. For this choice of 

C, a= UaVa 

— UaVa 

and 



b= Uav x — Uiv 3 . This vector is called the cross product . 

Now I know you're thinking: "There's no way in the world anyone 
can. remember that!!'' Aha--that's what you think!! Enter the 
world of determinants. I'll take you on a short detour, then come 
back to the cross product. 


DETERMINANTS 

A 2x2 determinant is four numbers arranged in the following 
pattern: 
la b| 


c d 


It is a number given by 


a b 

X 

c d 

* _ 'i 


= ad-bc 


A 3x3 determinant (and higher order determinants) are found in 
the fallowing manner: 



V w 


u w 


U V 

= a 


- b 


+ c 



y z 


X z 


x y 


x y z | 

a- fc-Q 
u v 

1 

x y | 

For the higher order determinants, follow the same pattern and 
continue to alternate signs as you go across the top row. 

Determinants have many interesting properties, and I have only 
given you the bare necessities to be able to calculate the value 
of a .determinant. But all the interesting stuff is not necessary 
for our immediate purposes. Now let's go back to the subject of 
cross products. 


®-b- -c 


u 

03 

jl V w 


1 

u y w 

Jc y z 


X ^ z 
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RETURN TO CROSS PRODUCTS 


Now recall, u = <u t ,us> ,u 3 > 

V = <Vi ,V 2 , Vs > 


We found a=u 2 v a - u 3 v 2 b=u 3 Vi - UiV 3 c=u t v 2 - u 2 vi. 




U2 

Ua 


Ui Ua 


U± Ua 

Using determinants: 

a = 



b= - 


C = 




Va 

Va 


V± Va 


Vi Va 


There is a nice "device" to get these 2x2 determinants by using 
the way we figure out 3x3 determinants. Toward this purpose, let 
me introduce a new notation for vectors: 

Note that any vector <p,q,r> can be written as 

p<l,0,0>+q<0,1,0>+q<0,0,1>. The vectors <1,0,0> and <0,1,0> and 
<0,0,1> are unit vectors pointing along the x-axis, y-axis, and z- 
axis respectively. These vectors have been given special names: . 

T = <1,0,0> 7 = <0,1,0> k = <0,0,1>. 


So our vector <p,q,r> = 
(Example: <3,4,5>=3i+2j+4k 

Ok, it now turns out that 


pi + qj + rk in our new notation. 


; <-2,1,3> = 

-2i+j+3k; etc.) 

- ' 

i ’ 

j 

k 

- 

U X V = 

Ux 

Ua 

Ua 

-- 


Vi 

Va 

Va 



To calculate, go across the top row as in our example in 
evaluating a 3x3 determinant. 



Ua 

Ua 


Ux Ua 


Ui Ua 



U X V = 

Va 

Va 

i - 

Vi Va 

j + 

Vi Va 

k 

Slick !!!! 


EXAMPLE: Let u = <2,1,3> and v = <l,4,-2> 



= (-2-12)i -(-4-3)j +(8-1)k 


= -14i +7j +7k = <-14,7,7> 
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PROPERTIES OF CROSS PRODUCTS 

Let me end by summarizing some properties of cross products. 
Some might be kind of hard to prove, but all should be fun to try 
to prove. 

(a) u x v = -(v x u) (ie. u x v is a vector in . opposite 
direction of v x u 

(b) u (v x w) = (u x v) w (note: there is no "normal" 

associative property of just cross products) 

(c) u x (v+w) = (u x v)+(u x w) (le. distributive prop.) 

(u+v) xw = (uxw) + (vxw) 

(d) (tu) x v = t(u x v)=u x (tv) 

(e) ju x v|= |u||v|sin0 

(f) u and v are parallel iff u x v = 0 

You now have most of the tools to tackle explanations and 
problems in 3-D space. Good luck! 
pr ob1ems 

1. Evaluate these determinants: 



3 2 


1 -2 


3 2 1 


1 2 1 

(a) 

1 4 

(b) 

4 7 

(c) 

1 4 1 

2 12 

.(d)- 

4 3 2 

2 11 


2. Exchange the 2nd and 3rd row in (d) above. Recalculate new 
determinant. Is the result "hinted at" in the properties above? 
Explain. 

3. Find u x v for these vectors: 

(a) u=<l,2,3> v=<l,-3,2> (b) u=<l / -l / 2> v=<-2,l,3> 

(c) u=<7,2,l> v=<-2,3,-2> 

4. Find i x j; j x k; i x k 
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chapter fourteen 
3—D TRANSFORMATIONS 


This follows along almost exactly the same as our previous two 
dimensional discussion. I'll only elaborate on any major 

differences. 

TRANSLATION 

To include translation as a matrix topic, we must go to 
homogeneous coordinates: (x,y,z)=(x,y,z,l) 

‘1 0 0 0 ! 


[x y z 1] 


0 10 0 

0 0 10 

H K L 1 


= [x+H y+K z+L 1] 


SCALING 


[x y z] 


Sx 0 0 

0 Sy 0 
0 0 S 3 


= [S*x S y y S*z] 


REFLECTION 
About x-y plane: 


[x y z] 


About x-z plane: 


About y-z plane: 

[x y z] 


10 0 
0 10 
0 0-1 

tx y z3 

-1 0 0 
0 10 
0 0 1 


[x y -z] 


10 0 
0-10 
0 0 1 


= [x -y z] 


[-x y z] 
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ROTATION 

We have 3 types of rotation-rotation about each of the axis. 




positioned so you can imagine looking down the x, z, and y axis 
respectively. Recall that positive rotation is counterclockwise 
from this view. Obviously (?) rotation around a particular axis 
doesn't affect that coordinate. (ie. rotation around the z-axis 
doesn't affect the z-coordinate of a point, etc.) 

The first two pictures above look exactly like the diagram - we 
worked from to derive the rotation formula in two dimension.- Thus 
we get: 


ROTATION ABOUT X-AXIS 

tx'y' z'] = [x y z] 


1 

0 

0 


0 0 
cos9 sin© 
-sin0 cos0 


X' = X 

y' = ycos0 - zsin0 
z' = ysin0 + zcos0 


ROTATION ABOUT Z-AXIS 


tx' y' z'] = [x y z] 


cos0 sin0 0 
-sin© cos0 0 
0 0 1 


x' = xcos© - ysin0 
y' = xsinS + ycos0 
z' = z 
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ROTATION ABOUT Y-AXIS 

This picture looks different!! We do get slightly different 
results. Suppose we start with a point (x,z) at some angle $, as 

Now rotate it by an angle 8 to get 

x=rsin$ and z=rcos4> 

x =rsin(<t>+0) and z =rcos(<I>+8) 

Simplifying x and z we get: 

x' = r(sin$cos8 + cos4>sin0) = rsin4>cose + rcos4>sin0 
= xcos0 + zsine 

z' = r(cos$cos0 - sin$sin8) = rcos$cos8- rsin$sin8 
= zcose - xsine 
= -xsine + zcos0 

So cos0 0 

[x' y'z']=[xyz] 0 1 

sine 0 

x' = xcose + zsine 

Y' = y - 

z' = -xsine + zcos8 

Ok, let's add these formulas for rotation to our previous 

ThreeD program in chapter 12. Rotation will be asked for in the 

input part of the program. 

The program is unsophisticated in that you can not change the 
order of the way the transformations are applied—and order does 
make a difference!! The program applies the transformations in 
the order that it asks for in the input! 

As always, feel free to improve this program. As with most 
learning, you learn best by doing! 


-sine 

0 

cose 


shown in the picture below. 
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A siaple deao of three-diaensional objects 
' With Perspective and Rotation added 
objectO: 

DATA 0,0,10,20,0,10,20,20,10,10,20,10,10,30,10,0,30,10,0,0,10,-1111,0,0 
DATA 0,0,0,20,0,0,20,20,0,10,20,0,10,30,0,0,30,0,0,0,0,-1111,0,0 
DATA 20,0,10,20,0,0,-1111,0,0 

DATA 10,20,10,10,20,0,-1111,0,0,10,30,10,10,30,0,-1111,0,0 
DATA 0,30,10,0,30,0,-1111,0,0,0,0,10,IT,0,0,-9999,0,0 
object 1: 

DATA 0,0,0,0,0,50,50,0,50,50,0,0,0,0,0,25,25,25,0,0,50,-1111,0,0 
DATA 50,0,0,25,25,25,50,0,50 
DATA -9999,0,0 

i 

LET xc=320 

LET yc=8Q ! Viewport [0,2*xc3x[0,2*ycl 
LET b=200 

LET d=100 ! window [-b,blx[-d,d] 

t 

DD 

CLS 

draw_axis 

set_up_textwindow 

INPUT “enter object # (0 or 1) to graph(-99 to end) »",object 

EXIT IF object=-99 
. ! 

set_up_textwindow 
LET~de=300 
LET ds=100 

PRINT "Scale of 1,no rotation, no translation, distance to eye is ";de 
PRINT "Distance to screen is "jds 
COLOR 1 

trig_calculations(0,0,0) 
plot_object(1,0,0,0) 

. fir 511iffle!=TRUE . ' 

' PAUSE 200 
• DO 

set_up_textwindow 

INPUT "Enter overall Scale factor—Sf »’,sf 

INPUT "Enter rotation angle about x-axis,y-axis,z-axis »",theta_x,theta_y,theta_z 
INPUT "Enter Translation in X and Y and Z direction »",h,k,l 
PRINT "Present distance to eye, distance to screen is "jde,ds 
INPUT "Enter distance to eye,distance to screen » ";de,ds 
IF firsttiae! THEN 
COLOR 0 ! Erase original figure 
trig_calculations(0,0,0) 
plot_object(1,0,0,0) 
draw axis 
firsttiae!=FALSE 
ENDIF 
COLOR 1 

trig_calculations(theta_x,theta_y,theta_z) 

plot_object(sf,h,k,l) 

set_up_textwindow 

PRINT "CHANGE OBJECTS OR CLEAR SCREEN?(y,or n) 
ansver$=INPUT$(l) 

EXIT IF UPPER$(answer$)=’Y" 

LOOP 

LOOP 
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PROCEDURE draw_axis 
COLOR 1 

DRAH 0,yc TO 2*xc,yc ! x axis 
DRAM xc,0 TO xc,2*yc ! y axis 
PRINT ATC77,11)jb 
PRINT AT(40,l);d 
RETURN 

PROCEDURE restore.data 
IF objects THEN 
RESTORE objectO 
ELSE 

RESTORE object 1 
ENDIF 
RETURN 

PROCEDURE set_up_textwindov 
FOR i=20 T0’25" 

PRINT AT(l,i)|SPACE$(73)j 
NEXT i 

PRINT AT(l,20)j 
RETURN 

PROCEDURE plot.object(scale,h,k,1) 
restore.data 
DO 

READ x,y,z 

scale.rotate.translate.object (scale, h,k,l) 
perspective.windoy 
IF NOT inside.figure! THEN 
PLOT x,y 
ENDIF 
DO 

READ x,y,z 

EXIT IF x=-9333 OR x=-llll 

scal e_r ot at e.t r ansi at e.ob jec t (scal e, h, k, 1) 

perspective.uindoy . 

IF NOT inside.figurs! THEN- 
DRAH TO- x,y ' ' 

ENDIF 

LOOP 

EXIT IF x=-3339 
LOOP 
RETURN 

PROCEDURE scale.rotate.translate.object(scale,h,k,1) 

LET x=x*scale 
LET y=y*scale 
LET z=z*scale 

’ rotation applied in order :abcut x-axis, then about y-axis 

' then about z-axis 

LET yy=y*ctx-z*stx 

LET z=y*stx+z*ctx 

LET y=yy 

LET xx=x*cty+z*sty 
LET z=-x*sty+z*cty 
LET x=xx 

LET xx=x*ctz-y*stz 
LET y=x*stz+y*ctz 
LET x=xx 

' Now translation 
LET x=x+h 
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LET y=y+k 
LET z=z+l 
RETURN 

PROCEDURE perspective_vindow 
LET inside figure!=FALSE 
LET dl=de-z 
IF dl>0 THEN 
LET x=x*(de-ds)/di 
LET y=y*(de-ds)/dl 
ELSE 

LET inside_figure!=TRUE 
ENDIF 

LET x=(x+b)/b*xc ! Nindowing function 
LET y=(d-y)/d*yc 
RETURN 

PROCEDURE trig_caIculations(ax,ay,az) 
stx=SIN(ax*PI/i30) 
ctx=C0S(ax*PI/180) 
sty=SIN(ay*PI/iSO) 
cty=C0S(ay*PI/100) 
stz=SIN(az*PI/180) 
ctz=C0S(az*PI/180) 

RETURN 
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AE>E>ENDI x 


ANSWERS 
PAGE 1-3 


1. PLOT 50,i 2. PLOT 110-1,50 3. PLOT i,10 

PLOT i, 100 
PLOT 100, i 
PLOT 10,i 


(|o,/o> 



PAGE 1-5 

DRAW 10,10 TO 10,100 TO 100,100 
FOR i=10 TO 90 STEP 10 
DRAW 10,i TO i+10,100 
NEXT i 




4. FOR m=l TO 100 STEP 2 
COLOR 1 

DRAW 10+m,10 TO 20+m,10 
DRAW TO 10+m,20 TO 10+m,10 
COLOR 0 

DRAW 10+m,10 TO 20+m,10 
- DRAW-TO 10+m,20 TO 10+m,10 
' NEXT m 

(I'm sure yours was more elegant!!) 


PAGE 2-3 

1. a) 7/4 b) -12/19 c) 11 d) 0 e) undefined 

2. (10,10),(15,12),(20,14)...many others 

3. (-1,0),(2,-5),(5,-10).. - 

4. They would have the same slopes.. 


PAGE 2-5 

1. a) y-6=-3(x-3) or y+3x=15 b) y-5=2/3(x+2) or 3y-2x=19 
c) y+4=-2/5(x+2) or 5y+2x=-24 

2. a) y=-3x+4 b) y=5 c) x=-3 

3. a) -3 b) 2/3 c) 0 d) 3 e) undefined f) -3/2 

4. a) y-4=-5/8(x-3) or 8y+5x=47 b) y-6=l/4(x-2) or 4y-x=22 
c) y+10=-8/7(x-3) or 7y+8x=-46 

5. y-4=-3/2(x-3) or 2y+3x=17 6. y-5=-5/3(x+4) or 3y+5x=-5 


PAGE 2-10 

1. a) <5,-4> b) <-9,12> c) <0,-3> d) <-6,14> e) <-3,-2> 

2. <-10,-6> 3. (7,-10) 4. (0,-8) 

5. 2y-x=5 and x-y=2—solve simultaneously and get x=9,y=7 
The points are (9,7) and (14,9). 

6 . Let u=<ul,u2>. u+u= <Ui , u 2 > + <ui, u 2 > = <ui+u 2/ Ui+Usa> = <2u * , 2u=> 

= 2<Ui, u=:> = 2u 

7. If t is positive, tu will be stretched(or shrunk) by a factor 
of t and tu will still be in the same direction. If t is 
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negative, tu will be in the opposite direction and stretched (or 

shrunk) by a factor of t . 

8. Any point 9. a 2 + b 5 * 


PAGE 2-12 

1. x=2+5t, y=-3+7tT 2. x=-3-47t, y=-7-t 

3. Horiz. line from (a,b) to (c,b) is x=a+(c-a)t, y=b 
Vert, line from (a,b) to (a,c) is x=a, y=b+(c-b)t 

4. If t=0 you get (xx,y»). If t=l you get (x a ,x = ) 

5. (7,4) and (15,1) 

6. If y=2, then 2=4-3t. Thus t=2/3. X=7+8(2/3)=37/3. 

The point is (37/3,2). If y=-2, the line segment does not intersect 

7. You get points on the "extended" line segment. 

8. Let (a,b) be the midpoint—for this point, t=l/2. 

So a=xi+(x 2 -xi) (l/2) = (xj.+X2t)/2. Similarly b=(y*+ya)/2 







1 

if 

b) - 3 


» c) 


W/,. 

- / 

-I 

1 

m 

3 


2. [0;639]x[0;199 ] . 

3. { (x,y,z): xe A and ye B and ze C } 


PAGE 3-3 

1. Suppose XL and YL are max. x & y values of screen. 
xv=(x+F)XL/2F yv=(G-y)YL/2G 

2. xv=same vv= v-C (YL-YS)+YS 

D-C 


PAGE - '4-3 

1. SIN(X), COS(X), SGN(X), EXP(X), etc. 

2. CHR$(X). 3. LEN(A$) 4. PLOT X,Y 

5. A subroutine often accepts certain numbers or strings(input) 
and performs a certain action(output). 

6 . a) 3/2 b) -1 c) 9 7. a) 12 b) 12 c) 4.44 

8 . a) (2,2) b) (4,4) c) (-3,-3) 

9. a) (6,1) b) (-6,7) c) (2p,4-p) 

10. a) (4,2) b) (-11,-3) c) (0.1,0.35) 

11 a)/T» 1.732 b) Undefined c) 4 a 2.588 

12. a) 0 b) * 0.9975 13. f(t)=(3-7t,4+t), Ostsl 

14. f (x) =2x+7 15. f (x,y)= /x 2 +y 2 

16. f( (xi ,y 2 ),(x 2 ,y a ) )=<x 2 -xi ,y 2 -y 2 > 


PAGE 5-2 

1. a) (x-5)* + (y-5)'* =100 b) (x+3) a +(y-6) a =4 

c) (x+4) 2 + (y+6) = =1/4 d) x a +y a =1 e) x 2 + y 2 

2. a) (0,2); r=6 b) (4,-3);r=*TT7 c ) (0,0 );r=/~8 

3. (x-5) a +(y-7)= =8 4. y=±3 x=±V 1.75 5. 


=25 

d) (-3.4); 
y=-3±v/T0' 


r = 8 


PAGE 5-8 

1. a) 3/5 b) 3/5 c) 4/5. d) 1 

appendix-2 



ST Hath & Graphics 


2. a) *0.8660 b) *0.9397 c) *0.6478 d) *-0.9528 e) *-0.8660 

3. 11.47 feet. 


PAGE 5-9 

1. a) *2.356 b) *1.221 c) *1.009 d) *2.124 e) *-1.356 

2. *11.05 3. a) 2sin(x)cos(x) b) cos a (x)-sin 2 (x) 

c) 2sin(x)cos(x) (there are many equiv. expressions) 
cos 2 (x)-sin^C x) 

4. cos^Cxj-sin^Cx) = cos = (x)-(l-cos =! (x) )= 2cos 2 (x)-l 

5. cosA=/3/2 ;V3/2 

6. sin(A+B)= 1/3 +/I5/6 ; cos(A+B)=/J/3 -/5/6 


1. 2y-3x=-ll 2. x=*+y 2 =l 


-T 




PAGE 6-5 

3. y=2x 3 -l 4. y 2 =4x 2 (1-x 2 ) 



1. a) x 5a +y 2 =2y b) x 2 +y 2 =2y+4x c) y=4x d) x 2 -y 2 =l 
e) x=2 f) y 2 =l+2x 

3. a) (2.6/1.5) b) (—2.2/—3.35) c) (1.6/1.2) d) (—2.85/1.45) 

4. a) (8.1/60.3*), (8.1/420.3*), (-8.1,240.3*) 

b) (8.9,116.6*), (8.9,-243.4*), (-8.9,296.6*) 

c) (8.1,240.3*), (8.1,-119.7*), (-8.1,60.3*) 

In each case, the initial angle was found from 
using the ATN(arctan or tan -1 ) function. 


tane =y/x -by 


1. A+B= fe 

3l ; AB= [*12 

13] 

r 0 “ o 

; BC= f8 11 

5 ’ 

5B= 

lo *5 ' 

; BA= 

7 

8 ' 

6 

9j Ll4 

21J 

1.26 34 

32 


15 25. 


14 

26 


b) No; Yes(BA=AB in part(a) ) 

c) C is a 2x3 while B is 2x3. The number of columns in A is not 
the same as the number of rows in B. 


2. 

'16 3 8 

AB= 

20 15 12 


.-17 -5 -5 

3. 

-3 4-2' 

-A= 

-5 -6 -2 


. 1 -2 3. 

4. a) 

AI=IA=A 

K interchanges 



9 

2 

0 * 


-13 

-32 -8 

BA= 

5 

6 

2 

A=>= 

43 

20 16 

L 

24 

-6 

11 


10 

10 11 


Define A-B = A+(-B) 


b) 

3 2 7 


3 7 2 

AK= 

2 14 

KA= 

2 7 2 


2 2 7. 


2 4 1 

2nd and 

; 3rd 

columns or 


3rd columns or 2nd and 3rd rows depending 


c) 

-3 7 2 

- c —u - 

-3 -7 -2 

AN= 

-2 4 1 

NA= 

2 4 1 


-2 7 2 


.2 7 2 

N changes the signs of the 1st col 


side you multiply. 
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1. 

2 0 

2. 

'.5 0* 

PAGE 

3 . ro i 

9-7 

; [. 707 -.707 

r ^ 

0 5 


0 .5. 

r 

L-i o. 

_ _-t 

_ [.707 .707 


1 

0 


0 

1 


10 -25 1 


-1 

0 

-20 


0 

1 

0 


-1 

0 

20 


0 

1 

0 


PAGE 11-5 

1. a) d=/x 2 +y a: +z i2 . b) /120 =2/30 c) /25 

3. a) (4.33,2.5,10) b) (4.95,-4.95,-2) c) (4.29,2.47,4.95) 
d) (-7.5,4.33,5) 

4. a) plane 5 units above the xy plane b) cone 

c) plane 4 units above xy plane d) vertical plane at 45* 
angle 


PAGE 13-3 

1. u*v =4; IuI= /74 ; IvI=/57 2. |vI= /84 

3. u*u =45 ; ( 1/YT5 )<5,4,2> 4. cos9= -49 _ 

/3T /S2 

5. 90 degrees 6. There are many—<1,1,1> will work.(u*v=0) 


PAGE 13-4 

1. z=0 2. -(x-2)+2(y+4)+7(z-3)=0 3. 2(x-3)-5(y-5)+3(z-7)=0 

4. -7(x-2)-(y-1)+10(z-3)=0 (or -7x-y+10z=15) 


PAGE 13-8 

1. a) 10 b) 15 c) 14 d) -1 

2. 1. Interchanging the rows give the opposite value. This is 
implied by property (a.) of cross products since if we calculate 
u x v and v x u, the vector is in the opposite direction. 

3. a) <13,1,-5> b) <-5,-7,-1> c) <-7,12,25> 

4. k ; i ; -j 
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append ± >c B 

Graphing program that was outlined in chapter 4. 

I'm sure you realize that your own efforts may not look anything 
like my program, and yet you may be perfectly correct. The main 
purpose of showing you MY program, is to give you some ideas in 
case you're stuck! 

' *** General graphing progras 
' *« Graphing 

' «* GFA BASIC VERSION 3 

' «* variables used 

’ xs,xl,ys,yl f a,b,c,d viewport and window 
' yi,v2,wl l w2 to speed calculations 
' ft function to be graphed 
’ answert response to input statement 
' x,y,xv,yv point in window and cooresp. point in viewport 

' HHHHHHHHHHHHHHHHHHHHHiiHHHHHHH 

LET f$="Cos(x)° i Edit for a different function * 

DEFFN f(x)=C0S(x) ! Edit for a different function * 

' **H*ti**m*ittmitii*m*M**f****t***m*t*t*i*miift 
DEFTEXT 2,16 ! Red, outlined text—“Bells and whistles* 

TEXT 32Q-LEN(f$)*4,8,UPPERt(f$) ! center title on screen 
’ default viewport and window 

DATA 0,639,0,159,-12,12,-2,2 - - ' - 

READ xs,xl,ys,yl,a,b,c,d . 

i I 

LET top=21 !top of text window 
REPEAT 
REPEAT 

eraseline(top) 

PRINT "The current viewport is C“jx5;","jxl;“lx[°;ysj",";ylj"]" 

eraselineftop+l) _ - . 

INPUT "Do you wish to change?(y/n) ",answer$ - , 

IF UPPER$(answer$)=“Y“ THEN _ " - ’/ 

eraselineftop+l) * 

INPUT "enter new viewport [Xs,Xl]xCYs,Y13 ",xs,xl,ys,yl 
eraselineftop+l) 

ENDIF 

UNTIL UPPER$(answer$)="N“ 

REPEAT 

eraseline(top+l) 

PRINT "The current window is [";aj","jb;"]x[";c;","jdj*l" 

INPUT "Do you wish to change?(y,n) "janswerf 
IF UPPERt(answer$)="Y" THEN 
eraseline(top+2) 

INPUT "Enter new window [A,Blx[C,D] ",a,b,c,d 
eraseline(top+2) 

ENDIF 

UNTIL UPPER$(answer$)="N" 

' following will speed up windowing function caluculations 

LET wl=b-a 

LET w2=d-c 

LET vl=xl-xs 

LET v2=ys-yl 

' Draw outline around viewport 
BOX xs,ys,xl,yl 

f 
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’ Draw in axis if they should be showing 
' Have adapted windowing function directly to give 
' the viewport coord, of axis 
IF a<0 AND b>0 THEN 
DRAH -a/wi*vl+xs,ys TO -a/wl*vl+xs,yl 
ENDIF 

IF c<0 AND d>0 THEN 
DRAH xs,-c/w2tv2+yl TO xl,-c/w2*v2+yl 
ENDIF 

' now draw curve 
FOR x=a TO b STEP wl/vl 
LET y=FN f(x) 

IF y<=d AND y>=c THEN Ipoint in window 
LET xv=(x-a)/wl*vl+xs 
LET yv=(y-c)/w2*v2+yl 
PLOT xv,yv 
ENDIF 
NETT x 

eraseline(top+2) 

INPUT "Repeat for another viewport/window?(y/n) ",answer$ 
eraseline(top+2) 

UNTIL UPPER$(answeri)="N" 

PROCEDURE eraseline(line) 

PRINT AT(l,line);SPACE$(73) 

PRINT ATCi,line); i reposition cursor at beginning of line 
RETURN 


Also included on the disk is a graphing program that 
the mouse to outline the viewport. It is called MOUSGRAF. 


uses 
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appendix C 

As promised in chapter 10, here is a general program to graph a 
function of two variables as a two-dimensional surface. The neat 
thing about this program is that it turned out to be a pretty 
short program. It doesn't have a lot of bells and whistles--no 
clipping, no true_ perspective, only graphs over the x-y plane in 
the first quadrant - . Yet it gives some nice results! 


The axis for this program are a bit 
unorthodox. I have the y-axis going back into 
the screen, as shown to the right. Use your 
right-hand rule to satisfy yourself that it is 
a right-handed coordinate system. The domain 
of a function of two variables can be thought 
of as the x-y plane. Of course we can't 
compute and plot points for all values in the 
x-y plane, so we lay out a checkerboard 
pattern and pick points at each intersection. 
The variables xdiv and ydiv specify how many 
subdivisions of the x and y axis we want. 



Below is a picture of the "window" and the viewport and the 
meaning of the various variables used in the program. 




Basically, we compute values of the points at the intersection 
of the checkerboard in the x-y plane, and save the values in a 
matrix Z. The points are then converted by the windowing function 
and plotted. The points are adjusted to the right and up to give 
the effect of depth and perspective. The fact that the y-axis 
looks like it is to the right and up, rather than being straight 
back into the screen is the basis of the adjustment. The 
variables dyxv and dyyv provide this adjustment. 
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' Deffn F(X,Y)=Sin(0.3*X)*Sin(0.3*Y) 

' Deffn F(X,Y)=(X/4) A 3-X/5*(Y/4) A 2 
DEFFN f(x,y)=COS(x*y)+l 

INPUT “Enter no. of div. in x-direction,y-direction “,xdiv,ydiv 

as 

DIN z(xdiv,ydiv) 

DATA Xb, 4 l Yo f 3,Za,4,Xc,10,Yc, 180,Xv,300,Yvx,300,Yvy, 100,Zv,50 
READ 

LET dx=xa/xdiv 
LET dy=ya/ydiv 
LET dxv=(xv-xc)/xdiv 
LET dyxv=(yvx-xc)/ydiv 
LET dyyv=(yvy-yc)/ydiv 
COLOR 3 

DRAH xc,yc TO xv,yc 
DRAW xc,yc TO yvx,yvy 

DRAH xc,yc TO xc,zv 
» 

FOR i=0 TO xdiv 
FOR j=0 TO ydiv 
z(i, j)=FN f(dx*i,dy*j) 

NEXT j 
NEXT i 

i 

COLOR 1 

FOR i=0 TO xdiv 
PLOT i*dxy,FN windowCz(i 1 0)) 

FOR j=i TO ydiv 

DRAN TO i*dxv+j*dyxv,FN windowCz(i,j))+j*dyyv 
NEXT j 
NEXT i 

FOR j=0 TO ydiv 

PLOT j*dyxv,FN windowCzC0 r jO)+j*dyyv_ .. ' • - 

FOR i=l TO xdiv .. . 

DRAH TO i*dxv+j*dyxv f FN windowCzCi.,j))+jtdyyv 
NEXT i 
NEXT j 

DEFFN window(zv)=zu/za*(zv-yc)+yc 
PRINT ATC1,25); 

INPUT “”,a$ 
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appendix E> 

You may wish to print out the listing o£ RPNSUB2.LST and 
.examine the routines. If added to your program, you may 
interactively enter in functions as an input string. If you 
remove the remark symbols from the first few lines of the 
"program" you can run a small demo program. There is also a 
graphing program, called INTERGRF that utilizes the subroutines. 

RPNSUB2.LST is a subroutine that converts a string that 
looks like a function into an RPN stack. The stack is then 
evaluated when values are desired. The program works pretty 
fast, not as fast as GFA's built-in machine language evaluation 
of functions, but we can't expect that. One of my favorite 
gripes is that none of the languages has a built-in interactive 
function input. I can't believe it'd be that difficult—the 
gut's of evaluating strings has to be already there in the 
language! Feel free to incorporate the subroutine in your own 
programs or improve them—I don't know how thoroughly the error 
checking will catch weird "mistake" entries. 
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